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1.  INTRODUCTION 

Prostate  cancer  (PCa)  is  one  of  the  most  common  cancers  for  males  in  the  United  States  and 
Canada,  with  245,000  estimated  cases  and  32,000  deaths  in  20151,2.  When  PCa  is  diagnosed  early, 
it  can  often  be  treated  by  surgery  or  radiation.  However,  in  about  30%  of  cases  when  PCa  recurs 
and/or  metastasizes,  the  disease  is  managed  primarily  by  the  hormone  therapy  that  reduces  the 
level  of  androgens  or  blocks  their  binding  to  the  androgen  receptor  (AR).  Unfortunately,  the 
effectiveness  of  this  therapy  is  temporary  due  to  emerging  resistance  mechanisms  related  to  AR. 
While  several  chemo-  and  immunotherapies  are  available  for  managing  such  advanced  PCa,  they 
offer  only  modest  survival  benefits.  Because  of  the  limited  effectiveness  and  significant  side 
effects  of  current  PCa  treatments,  there  is  a  pressing  demand  for  new  therapeutics  to  be  developed 
to  treat  aggressive  and  resistant  prostate  cancer.  With  the  accumulated  publicly  available  genomic 
data  and  recent  advancements  in  bioinformatic  and  cheminformatic  algorithms,  this  research 
project  aims  to  apply  comprehensive  analyses  of  cancer  genomes  and  transcriptomes  to  identify 
new  drug  targets,  followed  by  the  discovery  of  lead  compounds  through  virtual  screening 
techniques  and  biological  assays.  The  outcome  of  this  project  can  lead  to  new  generations  of  anti- 
PCa  drugs  to  treat  cases  where  the  cancer  cells  have  developed  resistance  to  existing  drugs,  by 
targeting  other  alternative  pathways. 
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2.  KEYWORDS 
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3.  ACCOMPLISHMENTS 


What  were  the  major  goals  of  the  project? 


Year  1  (August  1,  2014  to  July  31,  2015),  reporting  period 


Specific  Aim  1  -  Genomics:  identification  of  prostate 
cancer-related  aberrations  in  gene  regulation. 

Proposed  dates 
(month/year) 

Actual  completion 
dates  (month/year) 

Major  Goal  1:  Collect  available  genomic  data  from  prostate 
tumor  and  normal  cells. 

Aug/2014  - 
Nov/2014 

Aug/2014  - 
Nov/2014 

Milestone  1 :  Construction  of  a  shared  and  central  database  of 
genomic  data  for  prostate  cancer  research  at  Vancouver 
Prostate  Centre,  University  of  British  Columbia 

Nov/2014 

Nov/2014 

Major  Goal  2:  Analyze  the  genomic  data  to  identify  prostate 
cancer-specific  aberrations. 

Dec/2014  - 
May/2015 

Dec/2014  - 
May/2015 

Major  Goal  3:  Identify  new  drug  targets  in  prostate  cancer. 

June/2015  - 
July/2015 

June/2015  - 
July/2015 

Milestone  2:  Identification  of  three  drug  targets  in  prostate 

cancer. 

July/2015 

July/2015 

Milestone  3:  One  co-author  manuscript  on  bioinformatic 
analyses  of  genomic  data  and  drug  target  identification  in 
prostate  cancer. 

July/2015 

Dec/2014 

Year  2  (August  1,  2015  to  July  31,  2016) 

Specific  Aim  2  -  Cheminformatics:  development  of  new  Proposed  dates  Actual  completion 
classes  of  anti-prostate  cancer  drugs.  (month/year)  dates  (month/year) 

Major  Goal  4:  Construct  protein  structural  models  and  Aug/2015  - 

identify  potential  drug  binding  sites  for  each  of  three  targets.  Nov/2015 

Milestone  4:  Determine  binding  pockets  of  each  of  the  three  nov/2015 
targets  for  subsequent  virtual  screening. 


Major  Goal  5:  Virtual  compound  screening  against  the  Dec/2015  - 

binding  pocket  of  drug  target  #  1 .  Feb/2016 


Milestone  5:  Identify  200  compounds  with  high  binding 

potential  against  drug  target  #1  for  further  biological  Feb/2016 

validation. 


Major  Goal  6:  Virtual  compound  screening  against  the  Feb/2016  - 

binding  pocket  of  drug  target  #2.  Apr/20 1 6 
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Milestone  6:  Identify  200  compounds  with  high  binding 
potential  against  drug  target  #2  for  further  biological 
validation. 

Apr/2016 

Major  Goal  7:  Virtual  compound  screening  against  the 
binding  pocket  of  drug  target  #3. 

Apr/2016  - 
June/20 1 6 

Milestone  7:  Identify  200  compounds  with  high  binding 
potential  against  drug  target  #3  for  further  biological 
validation. 

June/2016 

Major  Goal  8:  Compound  validations  by  biological  assays. 

Aug/2015  - 
July/2016 

Milestone  8:  Identify  lead  compounds  with  high  biological 
activity  against  drug  target  #1,  #2  and  #3. 

June/20 1 6  - 
July/2016 

Milestone  9:  Up  to  three  co-author  manuscripts  on  lead 
compounds  against  drug  target  #1,  #2  and  #3. 

July/2016 

What  was  accomplished  under  these  goals? 

As  outlined  in  the  previous  section  and  in  the  approved  SOW,  the  three  major  goals  for  the  Specific 
Aim  1,  “Genomics:  identification  of  prostate  cancer-related  aberrations  in  gene  regulation”,  have 
been  accomplished  during  the  reporting  period,  Year  1  (August  1,  2014  to  July  31,  2015).  The 
three  accomplished  goals  are:  1)  collect  available  genomic  data  from  prostate  tumor  and  normal 
cells;  2)  analyze  the  genomic  data  to  identify  prostate  cancer-specific  aberrations;  and  3)  identify 
new  drug  targets  in  prostate  cancer. 

Major  Goal  1:  Collect  available  genomic  data  from  prostate  tumor  and  normal  cells. 

1)  Major  activities. 

Publicly  available  PCa  genomic  and  transcriptomic  datasets  from  high-throughput  experiments 
have  been  collected  from  a  number  of  sources  including  the  Gene  Expression  Omnibus  (GEO)3 
and  The  Cancer  Genome  Atlas  (TCGA)4.  The  public  available  datasets  were  downloaded  through 
the  corresponding  web  portals.  In  collaboration  with  Dr.  Anna  Lapuk,  in-house  high-quality 
datasets  previously  derived  from  the  patient  cohort  at  the  Vancouver  Prostate  Centre  (VPC), 
University  of  British  Columbia  (UBC)  have  been  collected.  Both  the  public  and  VPC  datasets  have 
been  stored  in  a  centralized  hard  drive  on  the  shared  VPC  network.  All  the  datasets  have  been 
checked  for  data  quality  and  consistency  (with  respect  to  their  original  publications).  These  data 
fdes  have  been  converted  and  stored  as  the  following  three  formats:  1)  tab-delimited  text  fdes  for 
easy  importing  to  other  analysis  programs;  2)  data  tables  in  the  MySQL  database  for  fast  data 
searching;  and  3)  data  objects  in  the  R  computational  environment  for  all  the  subsequent  statistical 
analyses.  A  number  of  scripts  based  on  the  PERL  programming  language  have  been  created  for 
data  processing  and  cross-referencing  among  the  difference  datasets. 
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2)  Specific  objectives 

The  main  objective  of  Major  Goal  1  is  to  create  a  comprehensive  collection  of  publicly  and 
internally  available  PCa  genomic  and  transcriptomic  datasets  not  only  for  the  analyses  performed 
in  this  project,  but  also  for  future  usages  by  other  colleagues  at  the  VPC,  UBC. 

3)  Significant  results/key  outcomes 

As  shown  in  Table  1,  this  comprehensive  collection  of  datasets  consists  of  2  DNA-Seq,  5  RNA- 
Seq,  12  microarrays  and  6  ChIP-Seq  datasets  from  a  total  of  1435  samples,  including  multiple 
prostate  cancer  cell  lines.  In  particular,  this  collection  also  contains  a  high-quality  in-house  dataset 
with  matched  DNA-Seq  and  RNA-Seq  data  from  26  prostate  tumours  and  5  normal  prostate 
samples  from  the  patients5. 


Table  1.  Prostate  cancer  (PCa)  genomic  and  transcriptomic  datasets  collect 

ted. 

Source/ Accession 

Data  type 

Total  number 
of  samples 

Sample  type 

PubMed  ID 

TCGA 

DNA-Seq, 
RNA-Seq,  SNP 
array 

550 

Patient  samples 

GEO/GSE38674 

DNA-Seq 

2 

LnCaP 

23245995 

VPC 

RNA-Seq 

31 

Patient  samples 

25155515 

GEO/GSE 14092 

RNA-Seq 

3 

LNCaP 

20478527 

RNA-Seq 

42 

Patient  samples 

22389870 

GEO/GSE34780 

RNA-Seq 

6 

LNCaP 

23034120 

GEO/GSE38676 

RNA-Seq 

2 

LNCaP 

23245995 

MSKCC/GEO/GSE2 1034 

microarray 

370 

Patient  samples 

20579941 

GEO/GSE60771 

microarray 

16 

VCaP 

25263440 

GEO/GSE58940 

microarray 

4 

VCaP 

GEO/GSE53994 

microarray 

15 

VCaP 

24591637 

GEO/GSE  15091 

microarray 

7 

PC3 

19668381 

GEO/GSE14028 

microarray 

15 

LNCaP,  VCaP 

20478527 

GEO/GSE28948, 

GSE35540 

microarray 

12 

VCaP 

22531786 

GEO/GSE36549 

microarray 

14 

LNCaP,  VCaP 

22710436 

GEO/GSE  15  805 

microarray 

155 

LNCaP 

22955617 

GEO/GSE28680 

microarray 

24 

Patient  samples 

23260764 

GEO/GSE39654 

microarray 

16 

LNCaP-lF5, 

VCaP 

23269278 
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GEO/GSE 11428 

microarray 

21 

LnCaP 

GEO/GSE 14092 

ChIP-Seq 

57 

LNCaP,  VCaP, 
patient  samples 

20478527 

GEO/GSE20042 

ChIP-Seq 

5 

LNCaP 

20208536 

GEO/GSE28950 

ChIP-Seq 

15 

VCaP 

22531786 

GEO/GSE38682, 
GSE38683,  GSE38684 

ChIP-Seq 

6 

LNCaP 

23245995 

GEO/GSE28219 

ChIP-Seq 

25 

LNCaP,  VCaP, 
patient  samples 

23260764 

GEO/GSE39879 

ChIP-Seq 

22 

LNCaP-lF5, 

VCaP 

23269278 

Major  Goal  2:  Analyze  the  genomic  data  to  identify  prostate  cancer-specific  aberrations. 

1)  Major  activities. 

In  order  to  identify  PCa-specific  aberrations,  the  datasets  collected  from  Major  Goal  1  have  been 
analyzed  by  using  the  following  computational  algorithms.  Single  nucleotide  variants  derived  from 
the  DNA-Seq  data  have  been  annotated  by  the  gSearch  program6,  and  copy  number  variations 
were  identified  by  the  use  of  BIC-seq7.  Gene  fusion  events  have  been  identified  by  using  DeFuse8 
on  the  RNA-Seq  data.  Differentially  expressed  genes  from  the  microarray  and  RNA-seq  data  were 
determined  by  using  t-test  for  two  group  comparison  (tumor  vs.  normal).  In  particular, 
overexpressed  genes  were  identified  from  each  dataset  by  the  following  steps:  1)  log2 
transformation;  2)  two  sample  t-test  between  tumour  and  normal  samples;  3)  multiple  testing 
correction  on  p-values9;  4)  selection  of  genes  with  corrected  (adjusted)  p-values  <  0.05;  and  5) 
among  those  with  significant  p-values,  selection  of  genes  with  fold-change  >  2  (tumor  vs.  normal). 
Transcription  factor  binding  data  from  ChIP-Seq  have  been  analyzed  by  using  the 
CompleteMOTIFs  pipeline10.  Additional  analyses  on  the  TCGA  dataset  were  performed  by  using 
the  computational  tools  available  at  cBioPortal11. 

2)  Specific  objectives 

The  main  objective  of  Major  Goal  2  is  to  perform  comprehensive  analyses  on  the  different  datasets 
(DNA-Seq,  RNA-Seq,  microarray,  ChIP-Seq)  collected  previously  from  Major  Goal  1.  The  results 
from  these  analyzes  have  been  used  to  identify  a  list  of  candidate  drug  targets  for  prostate  cancer 
(more  details  in  Major  Goal  3). 

3)  Significant  results/key  outcomes 

The  bioinformatic  analyses  performed  on  the  TCGA  prostate  cancer  dataset,  consisted  of  498  PCa 
patient  tumour  samples  and  52  normal  samples,  have  identified  a  list  of  frequent  genomic 
alterations,  as  shown  in  Table  2.  Among  the  top  alterations  are:  gene  fusions  involving  ERG, 
ETV 1  and  ETV4;  PTEN  deletion;  amplifications  of  genes  including  MY C,  AR,  and  PIK3CA.  The 
frequencies  of  these  genomic  alternations  identified  from  the  analyses  is  comparable  to  those 
reported  in  previous  studies12,13. 
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Table  2.  Frequent  genomic  alterations  in  PCa  patient  samples  (TCGA  dataset). 


Gene 

Type  of  alterations 

%  of  cases  in  the 
TCGA  dataset 

ERG 

Gene  fusion 

46% 

PTEN 

Deletion 

17% 

NKX3-1 

Deletion 

17% 

TMPRSS2 

Deletion 

16% 

ETV1 

Gene  fusion 

10% 

TP53 

Deletion 

8% 

MYC 

Amplification 

8% 

AR 

Amplification 

7% 

FOXA1 

Amplification 

6% 

ETV4 

Gene  fusion 

5% 

PIK3CA 

Amplification 

5% 

In  addition  to  the  genomic  alterations  listed  above,  gene  expression  analyses  performed  on  the 
different  datasets  collected  have  revealed  a  number  of  genes  that  are  overexpressed  in  PCa  patient 
samples.  Figure  1  illustrates  the  results  from  gene  expression  analyses  performed  on  the  two 
largest  datasets,  MSKCC14  and  TCGA4,  as  well  as  the  in-house  dataset  from  VPC5.  Using  a 
standard  but  stringent  bioinformatic  protocol,  a  total  of  169,  24,  744  genes  have  been  found  to 
overexpress  in  the  VPC,  MSKCC  and  TCGA  datasets,  respectively  (Figure  1A).  While  each 
dataset  features  a  different  list  of  genes,  ERG  is  the  only  common  gene  overexpressed  in  all  of  the 
three  datasets. 
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26  PCa  vs.  5  normal  samples  150  PCa  vs.  29  normal  samples 


To  identify  overexpressed  genes: 

1.  Log2  transformation 

2.  Two  sample  t-test  (tumour  vs.  normal) 

3.  Multiple  testing  correction 

4.  Adjusted  p-value  <  0.05 

5.  Fold  change  >  2 
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ERG  (VPC) 


"T 


’^T 

5  normal  26  tumour 


o 

•35  cnj 

to  T_ 
2? 

CL  O 


m 

C  CO 

0 

O) 

CN 

O) 

o 


ERG  (MSKCC) 


29  normal  150  tumour 


o 

to  cnj 
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0 
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ERG  (TCGA) 


52  normal  498  tumour 


Adjusted  p-value  =  2.62E-02 
Fold  change  =  3.29 


Adjusted  p-value  =  3.75E-15 
Fold  change  =  2.66 


Adjusted  p-value  =  4.53E-10 
Fold  change  =  2.90 


Figure  1.  ERG  is  one  of  the  most  commonly  overexpressed  genes  in  PCa  patient  samples. 

A)  A  Venn  diagram  that  shows  the  number  of  overexpressed  (i.e.  upregulated)  genes  from  each  of 
the  three  gene  expression  datasets:  Vancouver  Prostate  Centre  (VPC)5,  Memorial  Sloan-Kettering 
Cancer  Center  (MSKCC)14,  and  The  Cancer  Genome  Atlas  (TCGA)4,  based  on  a  bioinformatic 
protocol  as  listed  on  the  right  panel.  ERG  is  the  only  overexpressed  gene  from  intersecting  the 
three  datasets. 

B)  The  fold  changes  of  ERG  gene  expression  in  PCa  tumour  samples,  compared  to  normal 
samples,  range  from  2.66  to  3.29. 


To  further  explore  other  overexpressed  genes  (in  addition  to  ERG),  the  above  stringent 
bioinformatic  protocol  that  required  a  gene  to  be  overexpressed  in  all  three  datasets  has  been 
relaxed.  Figure  2  highlights  other  overexpressed  genes  in  the  TCGA  dataset,  which  represents  the 
largest  PCa  gene  expression  data  currently  available,  in  the  three  major  gene  families  of  interest: 
1)  ETS  factors,  2)  nuclear  receptors,  and  3)  kinases.  In  addition,  we  have  compared  the  gene 
expression  profiles  between  samples  with  low  ERG  expression  and  those  with  high  ERG 
expression  within  the  total  PCa  dataset  from  TCGA  (Figure  2A,  bottom  panel).  Notably,  ETV1 
and  ETV4  are  overexpressed  in  PCa  samples  where  ERG  is  expressed  at  lower  level.  This  gene 
overexpression  is  likely  due  to  the  gene  fusion  involving  ETV 1  and  ETV4,  and  it  suggests  that  the 
different  ETS  factors  (ERG,  ETV1  and  ETV4)  may  compensate  the  function  of  each  other. 
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A  ETS  family  B  Nuclear  Receptor  family  C  Kinase  family 


249  PCa  ERG  low  vs. 

249  PCa  ERG  high  ETS  protein  family 


*adj.  p-value  <  0.05,  fold-change  >=  1.3 
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Figure  2.  Overexpressed  genes  from  the  ETS,  Nuclear  Receptor  (NR)  and  Kinase  families  in 
PCa  patient  samples. 

A)  Top  panel:  A  list  of  ETS  genes  (ERG,  FEV,  SPDEF)  are  overexpressed  in  the  498  PCa  samples, 
compared  to  the  52  normal  samples  in  the  TCGA  dataset.  Bottom  panel:  Within  the  498  PCa 
samples,  4  ETS  genes  are  overexpressed  in  samples  that  have  lower  ERG  expression. 

B)  NR2E1  gene  from  the  nuclear  receptor  family  is  overexpressed  in  the  498  PCa  samples. 

C)  A  number  of  kinase  genes  are  overexpressed  in  the  498  PCa  smaples. 


To  gain  more  insights  on  DNA  binding  sites  by  the  ERG  transcription  factor  on  chromosomes  and 
nearby  genes,  several  publicly  available  ChIP-Seq  datasets  have  been  analyzed.  Figure  3  features 
the  result  from  analyzing  one  of  the  largest  ChIP-Seq  datasets15,  consisted  of  a  total  of  42,568 
DNA  sites  bound  by  the  ERG  transcription  factor  in  VCaP  cells  (a  PCa  cell  line  that  harbours  the 
ERG  gene  fusion  and  overexpresses  the  ERG  protein).  All  of  the  42,568  ERG-bound  sites  have 
been  analyzed  using  the  CompleteMOTIFs  program10,  and  Figure  3B  shows  the  top  nearby  genes 
ranked  by  the  number  of  ERG-bound  sites.  Interestingly,  the  gene  Androgen  Receptor  (AR)  is 
ranked  as  the  second.  Further  visualization  of  the  ERG-bound  sites  at  the  AR  locus,  using  the 
IJCSC  Genome  Browser16,  illustrated  that  the  sites,  with  various  amount  of  binding  intensity 
(ranged  from  0-250),  span  across  the  entire  AR  gene  including  both  upstream  and  downstream 
regions  (Figure  3A).  This  result  suggests  that  overproduction  of  the  ERG  protein  may  have  a 
significant  influence  on  the  AR  transcription  and  the  AR-regulated  pathways  in  prostate  cancer 
cells. 
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Entrez  gene  ID 

Gene  symbol 

Chromosome 

Transcription  start 

Transcription  end 

Number  of  ERG-bound  sites  (peaks) 

8660 

IRS2 

chrl3 

109204185 

109236915 

133 

367 

AR 

chrX 

66705408 

66860844 

121 

121793 

C13orfl6 

chrl3 

110771016 

110794595 

83 

60401 

EDA2R 

chrX 

65732207 

65752598 

57 

57507 

ZNF608 

chr5 

124000509 

124108704 

44 

6656 

SOX1 

chrl3 

111769914 

111774021 

40 

352887 

GYG2P1 

chrY 

13027923 

13043397 

38 

23769 

FLRT1 

chrll 

63627938 

63643221 

34 

283487 

C13orf29 

chrl3 

110314335 

110320656 

31 

55647 

RAB20 

chrl3 

109973414 

110012072 

29 

Figure  3.  Enrichment  of  ERG-bound  sites  at  the  Androgen  Receptor  (AR)  gene  locus  based 
on  ChIP-Seq  analyses. 

A)  Visualization  of  the  intensity  of  ERG-bound  sites  (y-axis)  at  the  AR  gene  locus  (x-axis)  on  the 
UCSC  Genome  Browser. 

B)  The  top  10  genes,  ranked  by  the  number  of  nearby  ERG-bound  sites. 


Major  Goal  3:  Identify  new  drug  targets  in  prostate  cancer. 

1)  Major  activities. 

By  analyzing  the  genomic  and  transcriptomic  datasets  as  summarized  in  Major  Goal  2,  we  have 
identified  a  number  of  target  genes  with  aberrations  such  as  genomic  alterations  and  gene 
overexpression  in  patient  PCa  samples.  In  order  to  perform  an  initial  assessment  of  structural 
availability  and  draggability  on  these  potential  targets,  the  protein  sequence  of  each  target  was 
queried  against  the  Protein  Data  Bank  (PDB)17  using  Protein  Blast18.  Protein  targets  with  resolved 
x-ray  crystal  or  NMR  structures  were  identified,  and  in  other  cases,  those  with  homologous 
structural  templates  (>  70%  sequence  identity)  suitable  for  subsequent  homology  modeling  were 
determined.  The  program,  Modeller  v9.1 1  (with  the  ‘allmodeT  protocol)19,  was  used  to  build  3D 
protein  models  for  the  potential  targets,  based  on  the  structural  templates  identified  from  PDB. 
Then  each  protein  model  was  subjected  to  the  Site  Finder  algorithm,  implemented  in  the  Molecular 
Operating  Environment  (MOE)  2012,  which  used  virtual  atomic  probes  to  search  the  protein 
surface  for  suitable  small  molecule  binding  pockets20.  Each  pocket  was  scored  based  on  its 
size/deepness,  solvent  exposure,  and  the  surrounding  amino  acid  residues. 

In  addition  to  the  two  criteria  as  described  above,  PCa-specific  aberrations  and  structural 
availability/draggability,  a  comprehensive  literature  has  been  performed  on  each  potential  target. 
In  particular,  we  examined  experimental  results  from  previous  studies  that  support  1)  gain  of 
cancerous  phenotypes  due  to  the  overexpression  of  the  target  and/or  2)  loss  of  cancerous  features 
due  to  the  knockdown  of  the  target.  All  these  information  were  presented  to  Dr.  Paul  Rennie 
(Mentor  of  this  project),  and  together  we  have  identified  a  total  of  three  candidate  drag  targets  that 
are  suitable  for  Specific  Aim  2  in  Year  2,  “Cheminformatics:  development  of  new  classes  of  anti¬ 
prostate  cancer  drags”. 
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2)  Specific  objectives 

The  main  objective  of  Major  Goal  3  is  to  narrow  the  list  of  potential  drug  targets  down  to  three 
candidates,  based  on  multiple  lines  of  supporting  evidence  based  on  1)  PCa-specific  aberrations  as 
identified  from  Major  Goal  2;  2)  existing  literature  supporting  the  oncogenic  functions  of  the  drug 
targets;  3)  expert  knowledge  from  Dr.  Rennie;  and  4)  structural  availability  and  draggability. 
Preliminary  protein  modeling  and  initial  binding  site  assessment  have  been  performed  on  each  of 
the  drag  target  candidates.  More  comprehensive  and  in-depth  structural  evaluations  based  on 
advanced  methods  such  as  molecular  dynamics  will  be  performed  in  Year  2  (Major  Goal  4: 
Construct  protein  structural  models  and  identify  potential  drag  binding  sites  for  each  of  three 
targets). 

3)  Significant  results/key  outcomes 

Initial  assessments  on  protein  structures  and  draggability  have  been  performed  on  a  number  of 
genes  with  PCa-specific  aberrations,  with  a  particular  emphasis  on  the  three  gene  families:  ETS 
factors,  nuclear  receptors  and  kinases.  Figure  4  shows  suitable  drag-binding  sites  identified  on  the 
surfaces  of  ERG,  ETV1  and  ETV4  protein  models.  To  further  examine  the  similarity  and 
difference  among  the  three  ETS  protein  sequences,  a  phylogenetic  tree  has  been  constructed  using 
a  multiple  sequence  alignment  program,  Clustal  Omega21.  As  shown  in  Figure  5,  the  ETS  family 
consists  of  a  total  of  28  members,  and  ERG  and  ETV1/4  are  clustered  to  the  ERG  and  PEA3 
subfamilies,  respectively.  In  addition,  the  sequence  alignment  indicates  that  ETV1/4  shares  61- 
65%  of  sequence  identity  with  ERG;  thus,  it  is  feasible  to  develop  specific  small  molecules 
targeting  each  of  the  three  ETS  factors  by  exploiting  their  sequence  and  structural  differences. 


ERG 


ETV1 


ETV4 


Figure  4.  Initial  assessment  of  druggable  sites  on  the  three  drug  target  candidates. 

Surfaces  of  predicted  protein  models  for  the  three  drag  target  candidates:  ERG,  ETV1  and  ETV4. 
Potential  drag-binding  sites  were  predicted  by  using  virtual  atomic  probes  (shown  as  red  spheres). 
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Figure  5.  ETS  phylogenetic  tree. 

A  phylogenetic  tree  of  the  28  members  in  the  ETS  family,  based  on  a  full-length  multiple  protein 
sequence  alignment. 


Combining  all  the  analysis  results  from  the  three  major  goals  (genomic  alteration,  gene 
overexpression,  protein  structure  availability  and  druggability,  as  well  as  supporting  evidence 
accumulated  in  the  literature12,13,22'25  (Table  3),  the  three  ETS  factors  (ERG,  ETVl  and  ETV4) 
have  been  identified  as  the  top  drug  target  candidates  to  be  moved  forward  into  Year  2  of  the 
project,  for  the  “Specific  Aim  2  -  Cheminformatics:  development  of  new  classes  of  anti-prostate 
cancer  drugs”. 


Table  3.  A  summary  of  the  three  drug  target  candidates  identil 

Red. 

Drug 

target 

Genomic 

alteration 

Gene 

overexpression 

Protein  structure 
availability 

Druggability 

Supporting 
evidence  in 
literature 

ERG 

Gene  fusion 

Fold-change  >  2 

yes 

yes 

yes 

ETVl 

Gene  fusion 

Fold-change  >1.3 

yes 

yes 

yes 

ETV4 

Gene  fusion 

Fold-change  >1.3 

yes 

yes 

yes 
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4)  Other  achievements 

The  bioinformatic  analyses  on  the  nuclear  receptors,  in  particular  the  orphan  nuclear  receptors,  as 
drug  targets  in  prostate  cancer  have  resulted  in  a  review  paper26.  A  manuscript  that  focuses  on  the 
ETS  factors  as  drug  targets  in  prostate  cancer  is  currently  in  preparation. 


What  opportunities  for  training  and  professional  development  has  the  project  provided? 

Training  activities: 

This  project  has  provided  me  (Dr.  Michael  Hsing,  postdoctoral  fellow)  a  great  amount  of 
opportunities  for  training  and  professional  development.  In  order  to  enhance  my  knowledge  in  the 
field  of  cancer  research  (and  in  particular,  prostate  cancer),  I  have  received  personalized  training 
from  my  mentor,  Dr.  Paul  Rennie,  who  has  over  30  years  of  research  experience  in  prostate  cancer 
and  more  than  180  peer-reviewed  publications  principally  in  this  field.  The  training  involved 
literature  review  of  previous  and  current  research  on  various  topics  in  prostate  cancer  and  other 
cancers  as  guided  by  Dr.  Rennie.  Cancer-related  journal  articles  and  textbooks  were  recommended 
by  Dr.  Rennie  for  studying.  My  knowledge  in  the  field  of  prostate  cancer  and  other  cancers  in 
general  have  been  examined  during  weekly  in-person  meetings  and  email  communication  with  Dr. 
Rennie.  Furthermore,  I  have  applied  this  knowledge  on  prostate  cancer  to  the  writing  of  a  review 
article26  in  this  field,  and  I  have  assisted  my  mentor,  Dr.  Paul  Rennie,  for  preparing  and  submitting 
two  research  grant  applications.  One  has  been  successfully  awarded  by  the  Pacific  Northwest 
Prostate  Cancer  SPORE  program,  while  the  other  has  been  submitted  to  the  Canadian  Cancer 
Society  Research  Institute.  In  addition,  I  have  gained  experience  on  analyzing  high-throughput 
sequencing  data  by  working  with  Dr.  Anna  Lapuk  and  other  bioinformatic  scientists  at  VPC,  UBC. 

Professional  development  activities: 

Vancouver  Prostate  Centre  (VPC),  UBC  is  located  at  the  largest  research  hub  in  Vancouver, 
adjacent  to  the  Vancouver  General  Hospital  and  its  affiliated  research  centers,  BC  Cancer  Agency 
(BCCA),  and  the  Genome  Science  Centre.  This  research  environment  greatly  facilitates  research 
communication  and  collaboration  among  not  only  local  scientists  but  also  international  researchers 
through  regular  seminars.  I  have  learned  about  the  latest  cancer  research  findings  and  interacted 
with  international  cancer  research  scientists  through  a  number  of  seminar  series  including  weekly 
Terry  Fox  seminars  hosted  at  VPC,  bi-weekly  seminars  at  BCCA,  and  monthly  SPORE  video 
conferences  hosted  by  the  Fred  Hutchinson  Cancer  Research  Center  at  Seattle,  Washington.  I  have 
presented  my  research  findings  at  regular  group  meetings  at  VPC,  and  I  attended  bi-weekly  journal 
clubs  at  VPC.  In  addition,  I  have  attended  and  presented  a  poster  at  an  international  conference, 
American  Association  for  Cancer  Research  (AACR)  Annual  Meeting  at  Philadelphia  during  April 
18-22,2014. 

How  were  the  results  disseminated  to  communities  of  interest? 

Nothing  to  Report. 

What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

In  accordance  with  the  approved  SOW,  I  have  accomplished  all  the  major  goals  and  tasks  in  this 
reporting  period  (Year  1:  Specific  Aim  1  -  Genomics:  identification  of  prostate  cancer-related 
aberrations  in  gene  regulation.).  A  total  of  three  drug  target  candidates  have  been  identified,  and 
the  project  is  on  track  for  the  next  reporting  period  (Year  2:  Specific  Aim  2  -  Cheminformatics: 
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development  of  new  classes  of  anti-prostate  cancer  drugs).  In  order  to  accomplish  the  goals  for 
Year  2, 1  will  follow  the  plan  as  described  in  the  approved  SOW.  In  particular,  the  following  goals 
will  be  achieved:  1)  construction  of  protein  structural  models  and  identification  of  potential  drug 
binding  sites  for  each  of  three  targets;  2)  virtual  compound  screening  against  the  binding  pockets 
for  each  of  the  three  drug  targets;  and  3)  compound  validations  by  biological  assays  (in 
collaboration  with  Dr.  LeBlanc). 
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4.  IMPACT 


What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

Prostate  cancer  (PCa)  is  one  of  the  most  common  cancers  for  males  in  the  United  States  and 
Canada,  with  245,000  estimated  cases  and  32,000  deaths  in  20151,2.  Major  issues  facing  PCa 
patients  are  cancer  recurrence  due  to  the  presence  of  distant  metastasis  at  time  of  surgery  and  the 
development  of  treatment  resistance.  The  proposed  project  is  aimed  to  address  this  urgent  issue  of 
metastasis  and  drug-resistance  in  prostate  cancer  by  using  the  combined  power  of  genomic 
analyses,  protein  structural  modeling,  computer-aided  drug  design  and  biological  screening  to 
develop  new  classes  of  anti-prostate  cancer  drugs,  targeting  previously  under-appreciated  proteins 
that  play  important  roles  in  cancer  progression. 

To  address  a  Prostate  Cancer  Research  Program  (PCRP)  overarching  challenge,  “Develop 
effective  treatments  and  address  mechanisms  of  resistance  for  men  with  high  risk  or  metastatic 
prostate  cancer”  with  a  focus  area,  “Therapy:  Identification  of  new  targets,  pathways,  and 
therapeutic  modalities”,  the  project  has  identified  three  drug  targets  based  on  large-scale  and 
comprehensive  analyses  of  available  genomic  and  transcriptomic  data  in  prostate  tumor  samples 
(as  a  result  of  the  current  reporting  period,  Year  1).  During  the  next  reporting  period  (Year  2),  each 
of  the  drug  targets  will  be  screened  rigorously  against  millions  of  virtual  chemical  models  to 
identify  candidate  molecules  with  the  best  structural  fit  in  the  binding  sites.  The  activity  of  these 
selected  compounds  will  then  be  validated  by  various  biological  assays  in  collaborations  with 
research  scientists  at  Vancouver  Prostate  Centre  (VPC),  University  of  British  Columbia. 

Although  it  is  a  lengthy  and  costly  process  to  take  a  drug  from  a  laboratory  to  a  patient  (with  an 
average  time  of  10-15  years  and  a  cost  up  to  billions  of  dollars),  the  time  and  cost  will  be 
significantly  reduced  by  the  computational  approaches  proposed  by  the  project  and  the  unique 
integration  at  VPC,  combining  drug  discovery  and  preclinical  development  under  a  single  research 
center.  We  anticipate  that  results  from  this  project  will  lead  to  new  drugs  that  can  be  used 
alternatively  or  synergistically  with  current  therapies  to  benefit  patients  with  the  most  deadly  and 
advanced  forms  of  prostate  cancer. 

What  was  the  impact  on  other  disciplines? 

Nothing  to  Report. 

What  was  the  impact  on  technology  transfer? 

Nothing  to  Report. 

What  was  the  impact  on  society  beyond  science  and  technology? 

Nothing  to  Report. 
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5.  CHANGES/PROBLEMS 


Changes  in  approach  and  reasons  for  change 

Nothing  to  Report. 

Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 

Nothing  to  Report. 

Changes  that  had  a  significant  impact  on  expenditures 

Nothing  to  Report. 

Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals,  biohazards, 
and/or  select  agents 

Nothing  to  Report. 

Significant  changes  in  use  or  care  of  human  subjects 

Nothing  to  Report. 

Significant  changes  in  use  or  care  of  vertebrate  animals. 

Nothing  to  Report. 

Significant  changes  in  use  of  biohazards  and/or  select  agents 

Nothing  to  Report. 
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6.  PRODUCTS 


Publications,  conference  papers,  and  presentations 
Journal  publications. 

Roshan-Moniri  M,  Hsing  M,  Butler  MS,  Cherkasov  A,  Rennie  PS.  Orphan  nuclear  receptors  as 
drug  targets  for  the  treatment  of  prostate  and  breast  cancers.  Cancer  Treat  Rev.  (2014).  40:  1 137- 
52.  (Published). 

Books  or  other  non-periodical,  one-time  publications. 

Nothing  to  Report. 

Other  publications,  conference  papers,  and  presentations. 

Oral  presentations: 

Hsing  M.  Development  of  small  molecule  inhibitors  of  ERG  to  treat  advanced  prostate  cancer. 
Prostate  Centre/  Terry  Fox  Seminar  Series.  Jim  Pattison  Pavilion,  Vancouver  General  Hospital, 
Vancouver,  Canada.  February  13,  2015.  (Local).  Acknowledgement  of  federal  support  (yes). 

Poster  presentations  and  abstracts: 

Roshan-Moniri  M,  Hsing  M,  Butler  MS,  Lau  D,  Axerio-Cilies  P,  Yen  P,  Kim  A,  Lien  S,  Mroczek 
M,  Ma  D,  Li  H,  Guo  Y,  Ban  F,  Ghaidi  F,  LeBlanc  E,  McIntosh  L,  Cox  ME,  Cherkasov  A,  and 
Rennie  P.  Therapeutic  targeting  of  ETS  factor  ERG  for  the  treatment  of  prostate  cancer.  AACR 
Annual  Meeting  2015.  Pennsylvania  Convention  Center,  Philadelphia,  USA.  April  18-22,  2015. 
(International).  Acknowledgement  of  federal  support  (yes). 

Butler  MS,  Hsing  M,  Roshan-Moniri  M,  Lau  D,  Yen  P,  Kim  A,  Lien  S,  Mroczek  M,  Ghaidi  F, 
LeBlanc  E,  McIntosh  L,  Cox  ME,  Cherkasov  A,  and  Rennie  PS.  Targeting  ETS  factor  ETV4  as  a 
novel  therapeutic  for  the  management  of  breast  and  prostate  cancer.  AACR  Annual  Meeting  2015. 
Pennsylvania  Convention  Center,  Philadelphia,  USA.  April  18-22,  2015.  (International). 
Acknowledgement  of  federal  support  (yes). 

Website(s)  or  other  Internet  site(s) 

Nothing  to  Report. 

Technologies  or  techniques 

Nothing  to  Report. 

Inventions,  patent  applications,  and/or  licenses 

Nothing  to  Report. 

Other  Products 

A  central  shared  database,  consisted  of  all  the  genomic  and  transcriptomic  datasets  collected 
during  this  reporting  period,  has  been  constructed  and  is  accessible  by  other  research  scientists  at 
the  Vancouver  Prostate  Centre,  University  of  British  Columbia. 
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7.  PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 
What  individuals  have  worked  on  the  project? 


Name: 

Michael  Hsing 

Project  Role: 

PI,  Postdoctoral  Fellow 

Researcher  Identifier 
(e.g.  ORCID  ID): 

Nearest  person  month 
worked: 

12 

Contribution  to  Project: 

Dr.  Hsing  has  performed  work  in  all  the  areas,  as  described  in  the 
Year  1  of  the  approved  SOW. 

Funding  Support: 

Name: 

Paul  Rennie 

Project  Role: 

Mentor,  Professor,  Director 

Researcher  Identifier 
(e.g.  ORCID  ID): 

Nearest  person  month 
worked: 

1 

Contribution  to  Project: 

Dr.  Rennie  has  provided  the  mentorship  and  guidance  to  Dr.  Hsing 
over  the  course  of  the  project. 

Funding  Support: 

Canadian  Cancer  Society  Research  Institute,  Canadian  Institute  of 
Health  Research,  Terry  Fox  Foundation 

Name: 

Anna  Lapuk 

Project  Role: 

Collaborator,  Research  Scientist 

Researcher  Identifier 
(e.g.  ORCID  ID): 

Nearest  person  month 
worked: 

1 

Contribution  to  Project: 

Dr.  Lapuk  has  provided  RNA-Seq  data  from  the  Vancouver 

Prostate  Centre,  University  of  British  Columbia. 

Funding  Support: 

Canadian  Cancer  Society  Research  Institute 

Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or  senior/key  personnel 
since  the  last  reporting  period? 

Nothing  to  Report. 
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What  other  organizations  were  involved  as  partners? 

Nothing  to  Report. 
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8.  SPECIAL  REPORTING  REQUIREMENTS 

None  required  (as  indicated  in  the  Award  document). 
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Nuclear  receptors  (NRs),  a  family  of  48  transcriptional  factors,  have  been  studied  intensively  for  their 
roles  in  cancer  development  and  progression.  The  presence  of  distinctive  ligand  binding  sites  capable 
of  interacting  with  small  molecules  has  made  NRs  attractive  targets  for  developing  cancer  therapeutics. 
In  particular,  a  number  of  drugs  have  been  developed  over  the  years  to  target  human  androgen-  and 
estrogen  receptors  for  the  treatment  of  prostate  cancer  and  breast  cancer.  In  contrast,  orphan  nuclear 
receptors  (ONRs),  which  in  many  cases  lack  known  biological  functions  or  ligands,  are  still  largely  under 
investigated.  This  review  is  a  summary  on  ONRs  that  have  been  implicated  in  prostate  and  breast  cancers, 
specifically  retinoic  acid-receptor- related  orphan  receptors  (RORs),  liver  X  receptors  (LXRs),  chicken  oval¬ 
bumin  upstream  promoter  transcription  factors  (COUP-TFs),  estrogen  related  receptors  (ERRs),  nerve 
growth  factor  IB-like  receptors,  and  “dosage-sensitive  sex  reversal,  adrenal  hypoplasia  critical  region, 
on  chromosome  X,  gene  1”  (DAX1).  Discovery  and  development  of  small  molecules  that  can  bind  at  var¬ 
ious  functional  sites  on  these  ONRs  will  help  determine  their  biological  functions.  In  addition,  these  mol¬ 
ecules  have  the  potential  to  act  as  prototypes  for  future  drug  development.  Ultimately,  the  therapeutic 
value  of  targeting  the  ONRs  may  go  well  beyond  prostate  and  breast  cancers. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


Introduction 

The  nuclear  receptor  (NR)  superfamily  is  composed  of  transcrip¬ 
tional  factors  that  enable  multicellular  organisms  to  regulate  gene 
expression  in  response  to  a  wide  variety  of  stimuli  from  develop¬ 
mental,  physiological  and  environmental  sources.  NR  actions  are 
modulated  by  several  endogenous  mechanisms  that  include:  (1) 
ligand  binding,  (2)  posttranslational  modification  of  amino  acid 
residues,  (3)  protein  dimerization,  (4)  nuclear  transfer,  (5)  pro¬ 
tein-protein  interactions  with  activators  and  repressors,  and  (6) 
cooperative  DNA  binding  with  other  transcription  factors  [1].  Of 
the  48  members  of  the  superfamily,  36  receptors  were  originally 
described  as  orphan  NRs  (ONRs)  because  their  exact  function  and 
native  ligands  were  undefined  [2].  As  a  result  of  recent  research 
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efforts,  21  of  these  ONRs  became  annotated  as  “adopted”  due  to 
the  discovery  of  natural  or  synthetic  ligands  (Fig.  1A). 

NRs  are  commonly  classified  into  7  subfamilies  based  on  their 
sequence  similarity  [3  and  have  been  studied  for  their  associations 
with  cancers  in  a  wide  range  of  human  organs  (Table  1).  However, 
there  is  a  disparity  of  research  efforts  dedicated  to  ONRs,  as  illus¬ 
trated  by  the  number  of  published  articles  in  PubMed  (Fig.  IB). 
For  example,  during  the  year  2013  alone,  543  and  1312  articles 
were  published  on  the  androgen  receptor  (AR)  and  estrogen  recep¬ 
tor  (ER)  in  the  context  of  their  corresponding  cancers.  By  compar¬ 
ison,  less  than  30  articles  were  published  in  the  same  year  on  the 
role  of  each  ONR  in  cancer,  as  highlighted  by  Fig.  1C.  This  imbal¬ 
ance  in  research  activity  is  also  mirrored  in  the  number  of  drugs 
and  chemicals  targeting  the  ONRs  (Fig.  1A). 

Conventionally,  targeting  transcriptional  factors  has  been  a 
challenging  task  in  drug  discovery  due  to  the  lack  of  well-defined 
active  sites  [4].  Thus,  the  unique  presence  of  the  ligand  binding 
domain  in  NRs  has  made  them  very  attractive  drug  targets,  and 
NRs  are  among  the  top  four  families  of  drug  targets  and  account 
for  approximately  13%  of  all  Food  and  Drug  Administration 
(FDA)-approved  drugs  [5].  To  gain  insight  into  the  number  of  exist¬ 
ing  small  molecules  targeting  the  NR  family,  each  NR  was  queried 
based  on  its  gene  symbol,  protein  name  and  protein  sequence 
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Fig.  1.  Comparison  of  existing  approved  drugs,  experimental  compounds  and  publications  for  NR  and  ONRs.  (A)  Each  NR  was  queried  for  the  number  of  approved  drugs  (the 
lefty-axis  shown  in  gray  bar)  in  DrugBank  [6],  and  for  the  number  of  published  compounds  (the  right y-axis  shown  in  black  bar)  with  experimentally-determined  binding 
affinities  in  BindingDB  [7]  in  October  2013.  The  NRs  are  arranged  according  to  the  subfamily  as  appeared  in  Table  1,  and  further  classified  as  orphan  or  adopted  ONR  according 
to  Shi  et  al.  [2]  Their  protein  sequences  were  compared  with  those  with  existing  experimental  structures  in  PDB  [25]  by  using  standard  protein  BLAST  [40].  NRs  that  have  no 
structural  template  with  sequence  identity  higher  than  70%  are  marked  by  **’.  The  number  of  articles  published  on  (B)  the  five  hormone  receptors  (ER,  AR,  PR,  GR  and  MR)  and 
PPARs  and  (C)  the  ONRs  discussed  in  this  review  paper,  were  extracted  from  PubMed  in  December  2013  using  a  search  query  consisted  of:  (Cancer  [Title/Abstract])  AND 
(“Name  of  the  nuclear  receptor” [Title/Abstract]).  Gene  names  and  synonyms  were  used  wherever  applicable. 


through  the  DrugBank  [6]  and  BindingDB  [7]  databases.  The 
DrugBank  contains  information  on  existing  drugs  and  their  direct 
targets,  including  1541  FDA-approved  small  molecule  therapeutics 
and  5082  experimental  drugs  that  interact  with  4323  proteins. 
BindingDB  currently  contains  over  one  million  binding  affinity 
records  corresponding  to  427,325  small  molecules  and  6589 
protein  targets,  extracted  from  the  literature.  As  illustrated  in 
Fig.  1,  five  hormone  NRs  (AR,  ER,  glucocorticoid  receptor  (GR), 


mineralocorticoid  receptor  (MR),  and  progesterone  receptor  (PR)) 
are  associated  with  a  large  number  of  synthetic  inhibitors  devel¬ 
oped  for  the  treatment  of  prostate-  and  breast  cancers  among  other 
diseases  [8,9].  Similarly,  some  adopted  ONR  such  as  peroxisome 
proliferator-activated  receptors  (PPARs)  have  received  consider¬ 
able  attention  and  have  clinically  approved  drugs  for  the  treatment 
of  diabetes  and  obesity  2].  In  contrast,  most  other  members  of  the 
ONR  family  have  fewer  associated  inhibitors  identified  due  to 
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Table  1 

Nuclear  receptors  in  cancers. 


Subfamily 

Group  name 

Abbreviation 

Cancer  association 

1.  Thyroid  hormone 

Thyroid  hormone  receptor 

TRot 

Liver  [166],  Kidney  [167],  Pituitary  [168],  Breast  [169],  Thyroid  [170] 

receptor-like 

(NR1A) 

TRp 

Retinoic  acid  receptor  (NR1B) 

RARoc 

Leukemia  [171],  Breast  [172],  Head  &Neck  [173],  Lung  [174],  Oral  tissue  [175],  Cervix  [176], 

RARp 

Ovary  [177],  Esophageal  [178],  Prostate  [179] 

RARy 

Peroxisome  proliferator- 

PPARoc 

Brain  [180],  Lung  [181],  Stomach  &  Colon  [182],  Liver  [183],  Pancreas  [184],  Bladder  [185], 

activated  receptor  (NR1C) 

PPAR-p/5 

PPARy 

Breast  [186],  Testis  [187],  Bone  [188] 

Rev-ErbA  (NR1D) 

Rev-ErbAoc 

Rev-ErbAp 

N/A 

RAR-related  orphan  receptor 

RORa 

Lymphoma  [189],  Breast  [73].  Prostate  [69] 

(NR1F) 

RORp 

RORy 

Liver  X  receptor-like  (NR1H) 

LXRot 

Prostate  [87],  Breast  [190],  Colon  [191],  Pancreas  [192],  Esophageal  [193],  Liver  [194] 

LXRp 

FXR 

Vitamin  D  receptor-like 

VDR 

Colon  [195],  Prostate  [196],  Breast  [197] 

(NR1I) 

PXR 

Colon  [198],  Breast  [199],  Ovary  [200] 

CAR 

Liver  [201] 

2.  Retinoid  X 

Hepatocyte  nuclear  factor 

HNF4a 

Liver  [202],  Colon  [203],  Breast  [204],  Prostate  [205] 

receptor-like 

(NR2A) 

HNF4y 

Retinoid  X  receptor  (NR2B) 

RXRot 

Breast  [146]  Colon  [206] 

RXRp 

RXRy 

Testicular  receptor  (NR2C) 

TR2 

N/A 

TR4 

TLX/PNR  (NR2E) 

TLX 

N/A 

PNR 

COUP/EAR  (NR2F) 

COUP-TFI 

COUP-TFII 

EAR-2 

Prostate  [106],  Lung  [207]  Breast  [108] 

3.  Estrogen 

Estrogen  receptor  (NR3A) 

ERa 

Breast  [9],  Cervix  [208],  Colon  [209],  Liver  [210],  Lung  [211],  Pancreas  [212],  Prostate  [213], 

receptor-like 

ERp 

Thyroid  [214],  Esophageal  [215],  Ovary  [216],  Adrenocortical  [217] 

Estrogen  related  receptor 

ERRot 

Prostate  [115],  Breast  [218] 

(NR3B) 

ERRp 

ERRy 

3-Ketosteroid  receptors 

GR 

Prostate  [219],  Breast  [220],  Lung  [221],  Pancreas  [222] 

(NR3C) 

MR 

N/A 

PR 

Ovary  [223],  Endometrial  [224],  Gallbladder  [225],  Breast  [226] 

AR 

Prostate  [8],  Breast  [227] 

4.  Nerve  growth 

NGFIB/NURR1/NOR1  (NR4A) 

NGFIB 

Breast  [228],  Colon  [229],  Prostate  [143] 

factor  IB-like 

NURR1 

NOR1 

5.  Steroidogenic 

SF1/LRH1  (NR5A) 

SF1 

Prostate  [230],  Adrenocortical  Tumors  [231] 

factor-like 

LRH-1 

Breast  [232],  Colon  [233] 

6.  Germ  cell  nuclear 

GCNF (NR6A) 

GCNF 

N/A 

factor-like 

0.  Miscellaneous 

DAX/SHP  (NROB) 

DAX1 

Bone  [234],  Breast  [155],  Prostate  [121] 

SHP 

Nomenclature  used  is  based  on  the  system  defined  by  the  Nuclear  Receptors  Nomenclature  Committee  [235]. 


difficulties  in  drug  development  and  toxicity,  as  well  as,  the  contin¬ 
ued  focus  on  the  identification  of  natural  ligands  [10].  For  example, 
the  role  of  retinoid  X  receptors  (RXRs)  as  heterodimeric  partners  of 
several  members  of  the  nuclear  receptor  family,  including  retinoic 
acid  receptor  (RAR),  makes  them  promising  drug  targets;  however, 
controversy  still  surrounds  whether  or  not  9-a’s-retinoic  acid  is  its 
true  endogenous  ligand  [11]. 

To  bring  attention  to  ONRs  as  prospective  novel  cancer  drug  tar¬ 
gets,  this  article  reviews  general  structural  features  and  biological 
functions  of  ONRs  and  discusses  some  of  the  most  promising  candi¬ 
dates  for  therapeutic  targeting  in  both  prostate  and  breast  cancers. 

Targeting  different  functional  sites  in  NRs  with  small  molecules 

All  NRs  are  modular  in  nature  and  composed  of  six  conserved 
regions,  designated  A-F,  which  make  up  the  following  domains: 
(A/B)  the  N-terminal  regulatory  domain  (NTD),  (C)  DNA  binding 


domain  (DBD),  (D)  hinge  region,  (E)  ligand  binding  domain  (LBD), 
and  (F)  C-terminal  domain.  A  phylogenetic  tree  was  constructed 
based  on  a  multiple  alignment  of  full-length  protein  sequences  of 
all  48  NRs  using  the  Clustal  Omega  [12]  software,  and  it  demon¬ 
strated  that  each  NR  can  be  further  grouped  into  the  corresponding 
subfamilies  based  on  sequence  similarities  (Fig.  2).  In  addition, 
Fig.  2  features  sequence  identity  between  each  NR  and  AR/ER  at 
the  DBD  and  LBD,  and  it  highlights  that  DBDs  are  more  conserved, 
sharing  ~50%  sequence  identity,  while  the  LBD  is  less  conserved 
(~25%  sequence  identity).  The  extent  of  sequence  conservation 
within  the  DBD  and  LBD  regions  are  further  depicted  as  a  heat 
map  in  Fig.  3A. 

Among  the  six  domains,  the  NTD  of  nuclear  receptors  contains  a 
transcriptional  activation  function-1  site  (AF-1),  and  its  activity 
does  not  directly  depend  on  the  presence  of  activating  ligands 
[13,14].  Although  the  transcriptional  activity  of  AF-1  is  limited,  it 
is  known  to  synergize  with  the  ligand-regulated  activation 
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SHP,  NR0B2  (none/none,  19/21%) 

— i 


{E 


—  TR -alpha,  NR1A1  (40/47%,  none/22'%) 

—  TR-beta,  NR1A2  (41/47%,  none/26%) 

-  RAR-alpha.  NR1&1  (43/57%,  24V23%) 

-  RAR-beta,  NR1B2  (43/57%,  24 */23%) 

-  RAR-gamma,  NR1B3  (43/60%,  24V26% ) 

—  HIMF4-alpha,  NR2A1  (48/47%,  26/26%) 

-  HNF4-gamma,  NR2A2  (50/49%,  25/27%) 
RXR-aipha,  NR2B1  (45/47%,  24/31%) 
RXR-gamma,  NR2B3  (46/49%,  24/32%) 

—  RXR-beta,  NR2&2  (45/49%,  25/32%) 

- «  ER-alpha,  NR3A1  (52/100%,  28/100%) 


ER-beta,  NR3A2  (58/84%,  30/60%) 
ERR^alpha,  NR3B1  (46/68%,  27/34%) 

I  ERR- beta,  NR3B2  (46770%,  28/37%) 


■m  ERR-gamma,  NR 3B3  (46/70%,  29/36%) 

-  GR,  NR3C1  (76/54%,  51/32%) 

- - MR,  NR3C2  (78/49%,  53/31%) 


—  PR,  NR3C3  (79/54%,  56/30%) 
-©  AR,  NR3C4  (100/52%,  100/28%) 


-V 


>  NGFIB,  NR4A1  (49/55%,  26/30%) 

I  NURR1,  NR4A2  (46/53%,  24/29%) 
«  NOR1,  NR4A3  (47/54%,  26/27%) 


-d 


»  COUP-TFI,  NR2F1  (45/49%,  24/26%) 

I  COUP-TFII.  NR2F2  (46/49%,  27/25%) 

EAR-2r  NR2F6  (49/52%,  3Q*/24%) 

- -  PPAR-alpha,  NR1C1  (42/49%,  27/31'%) 

- -  pp  A  R  -  beta! d  e  I  ta .  nr i  C2  (41/48%,  2t  *M6*%) 

- —  PPAR-gamma  NR1C3  (41/51%,  20/47*%) 

-  ReV-ErbA-alpha,  NRlDl  (38/49%,  25/23%) 

Rev-ErbA-betaL  NR1D2  (38/49%,  23/23%) 


- *  ROR-alpha,  NR1F1  (43/48%,  29*/23%) 

ROR-beta.  NR1F2  (43/47%.  28'/26*%) 

- •  ROR-garnma,  NR1F3  (43/44%,  3 V/24*%) 

LXR-alpha,  NR1H3  (46/52%,  25/21 U) 

LX  R- beta,  NR1H2  (48/54%,  22/20%) 

-  FXR,  NR1H4  (48/51%,  28/23%) 


^ - £ 


—  VDRf  NRltl  (43/47%,  27*1 36*%) 

—  PXR,  NR112  (45/48%,  21/31*%) 

—  CAR  t  NR1 13  ( 44/45%.  28*/none%) 
SF1,  NR5A1  (49/50%,  28/28%) 

—  LRH-1,  NR5A2  (47/51%,  23/28%) 

—  TR2,  NR2C1  (50/47%,  24/25%) 
TR4,  NR2C2  (51/49%,  2SV22%) 


-  GCNF,  NR6A1  (45/49%,  24/25%) 

—  PNR,  NR2E3  (43/44%,  31/27%) 

-  TLXr  NR2E1  (44/51%,  25/31%) 

\  DAX1.  NR0B1  (17Vnone%,  none/22*%) 


Fig.  2.  A  phylogenetic  tree  of  NRs.  A  phylogenetic  tree  was  constructed  based  on  a  multiple  alignment  of  full-length  protein  sequences  of  all  48  NRs  using  Clustal  Omega  [12]. 
The  ONRs  featured  in  the  review  paper  are  highlighted  in  red,  and  the  two  nuclear  receptors,  AR  and  ERa,  are  indicated  in  green  as  the  reference.  The  sequence  identity  of  each 
NR  is  compared  to  AR  and  ERot  at  the  DBD  and  LBD,  with  the  domain  boundary  positions  defined  by  NCBI  Conserved  Domain  Database  (CDD)  [161].  The  following  format  is 
used  (a/b%,  c/d%):  a  =  %  identity  vs.  AR-DBD,  b  =  %  identity  vs  ERoc-DBD,  c  =  %  identity  vs  AR-LBD,  and  d  =  %  identity  vs  ERoc-LBD.  *  =  query  alignment  coverage  is  less  than  50%. 


function-2  (AF-2)  regulatory  site  located  within  the  LBD  to  regulate 
gene  expression.  Significant  variability  in  the  NTD  can  be  observed 
even  within  a  single  type  of  NR  where  different  sequences  can  be 
derived  from  alternative  splicing,  promoters,  and  transcription 
start  sites.  One  of  the  best  examples  of  such  variability  is  found 
with  RAR  where  three  gene  templates,  designated  oe,  p  and  y  yield 
eight  alternative  protein  variants  [15,16]. 

The  DBD  has  the  highest  sequence  conservation  among  the  NRs, 
and  it  mediates  receptor  binding  to  response  elements  (REs) 
located  in  the  enhancers  and  promoters  upstream  of  NR-regulated 
genes.  The  receptors  can  bind  to  REs  as  a  monomer,  homodimer  or 
heterodimer,  and  the  two  zinc  finger  motifs  within  the  DBD  aid  RE 
binding  and  receptor  dimerization  [17].  All  NRs  contain  a  DBD, 
with  the  exception  of  dosage-sensitive  sex  reversal,  adrenal  hypo¬ 
plasia  critical  region,  on  chromosome  X,  gene  1  (DAX1)  and  small 
heterodimer  partner  (SHP),  both  of  which  lack  a  DNA  binding 
domain  [18,19].  The  absence  of  DBD  on  DAX1  and  SHP  can  also 
be  demonstrated  by  the  lack  of  sequence  alignment  in  Fig.  3A. 

The  hinge  region  is  the  connecting  bridge  between  the  DBD  and 
LBD  sections.  The  sequences  of  the  hinge  region  vary  significantly 


among  NRs  and  contribute  to  flexible  rotations  of  the  DBD, 
enabling  direct  or  inverted  binding  and  dimerization  of  NRs  onto 
the  HREs  [20].  There  is  some  evidence  that  this  region  may  also 
serve  as  a  binding  site  for  co-repressor  proteins  [21,22]. 

The  LBD  has  only  moderate  sequence  conservation  among  NRs, 
but  possess  highly  similar  three-dimensional  organizations.  The 
LBD  is  composed  of  11-12  alpha  helices,  which  can  fold  into  a  com¬ 
pact  formation  with  parallel  and  anti-parallel  orientation  (Fig.  3B). 
This  domain  is  associated  with  multiple  functions  including  ligand 
attachment  through  the  ligand-binding  site  (LBS),  interactions 
with  other  proteins  at  the  AF-2  and  binding  function-3  (BF-3)  co¬ 
regulator  binding  sites,  protein  dimerization,  and  nuclear  translo¬ 
cation  [23].  For  more  in-depth  discussions  on  AF-2  and  BF-3  sites 
in  NRs,  we  refer  readers  to  the  two  recent  review  papers  by  Buzon 
et  al.  [24]  and  Caboni  and  Lloyd  [17]. 

In  additional  to  high-throughput  screening  and  other  conven¬ 
tional  approaches,  computer-aided  drug  design,  guided  by  avail¬ 
able  protein  structures  in  the  Protein  Data  Bank  (PDB)  [25]  and 
synthetic  compounds  from  databases  such  as  ZINC  [26]  can 
provide  a  cost-  and  time-effective  strategy  to  search  for  new  drug 
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Fig.  3.  Sequence  conservation  and  a  general  3D  structural  model  of  protein  domains  in  nuclear  receptors.  (A)  A  heat  map  at  the  top  panel  illustrates  a  multiple  protein 
sequence  alignment  among  the  48  NRs  at  the  DNA  binding  domain  (DBD)  and  ligand  binding  domain  (LBD),  generated  by  using  Clustal  Omega  [12]  and  visualized  in  Jalview 
[162].  The  domain  boundary  positions  are  defined  according  to  NCBI  CDD  [161  ].  A  darker  color  represents  a  residue  that  matches  the  consensus  sequence  (more  conserved)  at 
a  position,  while  a  lighter  color  represents  a  different  residue  but  with  a  positive  Blosum62  score  (less  conserved).  A  non-conserved  residue  or  a  gap  is  colored  in  white.  The 
48  NRs  are  arranged  from  top  to  bottom  (y-axis)  in  the  same  order  as  they  appear  in  Table  1.  The  aligned  residues  are  arranged  from  the  N-terminus  to  the  C-terminus  on  the 
x-axis.  Residues  within  the  LBS,  AF-2  and  BF-3  sites  are  highlighted  as  black  bars  on  top  of  the  LBD  heat  map,  using  the  AR  sequence  as  the  reference  positions  defined  by  NCBI 
CDD  [161]  and  Buzon  et  al.  [24].  (B)  The  bottom  panel  shows  a  general  3D  structural  model,  reconstructed  by  combining  several  crystal  structures  of  individual  domains, 
which  include:  (1)  DBD  (in  pink  ribbon)  bound  at  DNA  (PDB  ID:  1R4I)  [163];  (2)  the  relative  orientation  among  DBD,  the  hinge  region  and  LBD  (PDB  ID:  4IQR)  [164];  (3)  LBD 
(in  orange  ribbon)  and  the  ligand  testosterone  bound  at  the  ligand-binding  site  (PDB  ID:  2AM9)  [165];  (4)  LBD  with  a  small  molecule  inhibitor  bound  at  the  AF-2  site  (PDB  ID: 
2YHD)  [35];  and  (5)  LBD  with  a  small  molecule  inhibitor  bound  at  the  BF-3  site  (PDB  ID:  4HLW)  [37].  No  structural  model  of  the  N-terminal  regions  is  available.  Protein 
models  were  visualized  by  using  the  Molecular  Operating  Environment  (MOE). 


candidates  and  has  been  successfully  applied  to  a  number  of  NRs 
[17,23].  For  instance,  virtual  screening  methods  have  been  used 
to  discover  agonists  and  antagonists  targeting  the  ligand  binding 
domains  of  PPAR  [27],  RAR  [28],  estrogen  related  receptor  (ERR) 
[29  and  Liver  receptor  homolog  1  (LRH-1)  [30].  While  many  exist¬ 
ing  drugs  targeted  the  ligand  binding  site  within  the  LBD  of  NRs, 
such  as  Enzalutamide  [31]  and  Bicalutamide  [32]  for  AR  and 
Tamoxifen  33]  and  Fulvestrant  [34]  for  ER,  there  are  several  new 
inhibitors  that  have  been  developed  to  target  protein-protein 
interactions  at  alternative  sites,  as  exemplified  by  AF-2  [35]  and 
BF-3  [36,37]  of  the  androgen  receptor  (Table  3).  These  alternative 
sites  (depicted  in  Fig.  3B)  are  required  for  physical  interactions 


with  protein  co-factors;  thus,  binding  of  small  molecules  at  such 
sites  can  block  co-factor  interactions  and  provide  additional  means 
for  inhibiting  NRs,  particularly  in  cases  where  drug  resistance  has 
occurred  due  to  mutations  in  the  LBS.  Furthermore,  it  has  been 
shown  recently  as  a  proof  of  concept  that  the  protein-DNA  inter¬ 
face  on  the  DBD  of  AR  can  be  targeted  by  small  molecules 
[38,39],  and  this  work  opens  up  a  new  strategy  to  inhibit  different 
protein  variants  that  may  differ  in  LBD  but  share  the  same  DBD 
within  the  AR  or  in  other  NR  subfamilies. 

To  evaluate  the  “druggability”  and  the  number  of  available 
structural  templates  of  each  NR,  protein  sequences  were  aligned 
with  those  with  existing  protein  structures  in  PDB  [25]  by  using 
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Table  2 

A  list  of  orphan  nuclear  receptors  featured  in  this  review  paper. 


Name 

UniProt 

ID 

Endogenous  ligands 

Interactors 

(BioGRID) 

Protein  structures  (PDB  ID) 

LBS/AF-2/BF-3  vs.  AR 
(sequence  identity, 
similarity) 

RORoc 

P35398 

Cholesterol,  cholesterol  sulfate, 
hydroxycholesterol,  ketocholesterol, 
epoxycholesterol,  ruscogenin  [64] 

24  Interactors 
(ID:  112022) 

1N83,  1S0X 

6/27/17%,  24/64/25% 

ROR[3 

Q92753 

Stearic  acid,  retinoic  acid  [64] 

14  Interactors 
(ID:  112023) 

1K4W*,  1NQ7*,  1N4H* 

6/27/17%,  24/64/25% 

RORy 

P51449 

Hydroxycholesterol,  ketocholesterol, 
epoxycholesterol  [64] 

3  Interactors 
(ID:  112024) 

3 BOW,  3KYT,  3L0J,  3L0L,  4NB6,  4NIE 

12/27/25%,  29/64/33% 

LXRcx 

Q.1 3 1 33 

Oxysterols  [90] 

37  Interactors 
(ID:  115373) 

1UHL,  3IPQ,  3IPS,  3IPU 

18/36/25%,  35/64/25% 

LXR[3 

P55055 

Oxysterols  [90] 

54  Interactors 
(ID:  113222) 

1P8D,  1PQ6,  1PQ9,  1PQC,  1UPV,  1UPW,  3KFC, 
3L0E,  4DK7,  4DK8 

18/36/25%,  35/64/25% 

COUP- 

TFI 

P10589 

Orphan 

20  Interactors 
(ID:  112883) 

2EBL 

18/9/17%,  35/73/25% 

COUP- 

TFII 

P24468 

Orphan 

30  Interactors 
(ID:  112884) 

3CJW 

18/9/17%,  35/73/25% 

ERRoc 

PI 1474 

Orphan 

44  Interactors 
(ID:  108405) 

1XB7,  2PJL,  3D24,  3K6P 

35/45/25%,  76/91/33% 

ERR[3 

095718 

Orphan 

7  Interactors 
(ID:  108406) 

1L01 

35/45/25%,  71/91/42% 

ERRy 

P62508 

Orphan 

24  Interactors 
(ID:  108407) 

1KV6,  1TFC,  1VJB,  2E2R,  2EWP,  2GP7,  2GPO, 
2GPP,  2GPU,  2GPV,  2P7A,  2P7G,  2P7Z,  2ZAS, 
2ZBS,  2ZKC 

35/45/25%,  71/91/33% 

NGFIB 

P22736 

Orphan 

84  Interactors 
(ID:  109407) 

2QW4,  3V3E,  3V3Q,  4JGV,  4KZI,  4KZJ,  4KZM 

12/9/33%,  47/45/33% 

NURR1 

P43354 

Orphan 

14  Interactors 
(ID:  110983) 

10VL 

18/0/33%,  47/45/42% 

NOR1 

092570 

Orphan 

4  Interactors 
(ID:  113713) 

1CIT* 

18/0/33%,  47/45/42% 

DAX1 

P51843 

Orphan 

15  Interactors 
(ID:  106695) 

3F5C** 

6/27/17%,  41/100/33% 

The  lack  of  endogenous  ligands  in  ONRs  was  confirmed  by  querying  the  IUPHAR/BPS  database  [236]  in  October  2014.  BioGRID:  the  number  of  interactors  queried  in  October 
2014.  PDB:  available  protein  structures  queried  in  July  2014.  Protein  structure:  ^homologous  protein  structures  with  ^  70%  sequence  identity.  **homologous  protein 
structures  with  <70%  sequence  identity.  Sequence  similarity  was  calculated  based  on  the  Blosum62  matrix  vs.  AR  sequences.  Residue  positions  within  each  functional  site, 
LBS  (n  =  17),  AF-2  (n  =  11)  and  BF-3  (n  =  12),  were  defined  based  on  the  AR  protein,  according  to  NCBI  Conserved  Domain  Database  (CDD)  [161]  and  Buzon  et  al.  [24]. 


a  standard  BLAST  algorithm  [40].  As  highlighted  in  Fig.  1  and 
Table  2,  the  majority  of  NRs  and  ONRs  have  protein  structural  tem¬ 
plates  with  sequence  identity  more  than  70%,  which  is  sufficient 
for  homology  modeling  as  a  basis  for  virtual  screening.  The 
multiple  sequence  alignment  as  represented  as  a  heat  map  in 
Fig.  3A  illustrates  that  the  positions  of  amino  acids  at  LBS,  AF-2 
and  BF-3  sites  overlap  with  highly  conserved  regions  within  LBD. 
This  is  in  accordance  with  previous  reports  [24]  indicating  that 
these  functional  sites  are  often  conserved  among  many  NRs  includ¬ 
ing  members  of  the  ONR  family  [41  -45  ]  and  as  such  may  be  target- 
able.  To  facilitate  future  drug  design  targeting  the  ONRs,  Table  2 
shows  a  list  of  ONRs  featured  in  this  review  paper  and  their  anno¬ 
tations  including  accession  numbers  in  the  UniProt  database  [46], 
endogenous  ligands,  the  number  of  interactors  listed  in  the 
BioGRID  database  [47],  and  available  protein  structures  in  PDB. 
To  determine  the  extent  of  sequence  conservation  at  the  LBS  and 
possible  presence  of  AF-2  and  BF-3,  the  corresponding  residues 
of  each  ONR  from  the  multiple  sequence  alignment  were  compared 
to  those  of  the  androgen  receptor,  and  the  sequence  identity/ 
similarity  values  were  listed  in  Table  2.  Interestingly,  the  AF-2  site 
of  RORs,  LXRs,  COUP-TFs,  ERRs  and  DAX1  share  moderate  to  high 
sequence  similarities  with  the  AR  (64%,  64%,  73%,  91%  and  100% 
respectively),  and  previous  studies  have  confirmed  the  presence 
of  AF-2  in  RORs  and  ERRs  [48,49]. 

In  the  next  sections,  we  focus  on  these  ONRs  as  potential  drug 
targets  in  prostate  and  breast  cancers.  For  each  ONR,  we  review 
its  functions  in  normal  physiological  and  cancer-related  condi¬ 
tions.  The  druggability  of  each  highlighted  ONR  and  known  small 
molecule  ligands  are  also  discussed.  Table  3  summarizes  a  selected 
list  of  existing  drugs  and  experimental  compounds  targeting  each 


ONR,  while  synthetic  compounds  targeting  AR  and  ER  are  also 
included  for  reference.  For  the  roles  of  other  NRs  and  ONRs  in  dif¬ 
ferent  cancers,  we  refer  readers  to  a  list  of  representative  articles 
summarized  in  Table  1. 

Targeting  ONRs  in  prostate  and  breast  cancers 

Prostate  cancer  is  one  of  the  most  common  cancers  in  men,  with 
an  estimated  238,590  newly  diagnosed  cases  and  29,720  deaths  in 
2013  in  the  United  States  alone  [50].  As  a  result  of  improved 
screening  procedures,  prostate  cancers  are  currently  being  diag¬ 
nosed  earlier,  when  the  disease  is  localized  to  the  prostate;  if 
detected  early,  the  treatment  options  of  radical  prostatectomy 
and/or  radiation  therapy  are  potentially  curative.  However,  the 
main  treatment  options  for  advanced  or  metastatic  prostate  can¬ 
cers  are  hormonal  therapies. 

Prostate  cancer,  like  the  non-malignant  prostate,  is  dependent 
on  androgens  for  proliferation  and  survival.  Thus,  hormonal  thera¬ 
pies  that  aim  to  either  reduce  androgen  levels  (via  surgical  or 
chemical  castration)  or  inhibit  the  AR  protein  (by  small  molecules) 
are  used  to  slow  tumor  growth.  Unfortunately,  despite  initial 
response,  such  treatments  will  eventually  fail  and  the  cancer 
recurs  in  an  incurable,  castration-resistant  form  [51,52].  This  treat¬ 
ment  resistance  may  be  a  result  of  several  mechanisms  including 
amplification  of  AR  expression,  mutations  at  the  androgen-binding 
site  of  the  AR  and  elevated  production  of  AR  splice  variants  [53]. 
Thus,  the  development  of  new  therapeutics  to  combat  metastatic 
castration-resistant  prostate  cancer  remains  a  major  task. 

Breast  cancer  is  the  leading  cause  of  cancer-related  death  in 
women  worldwide  [54].  In  2013  alone,  there  were  an  estimated 
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Table  3 

A  selected  list  of  compounds  targeting  nuclear  receptors. 


Compound  name 


Chemical  structure 


Annotation  (protein  target/accession  number/status/activity) 


Enzalutamide 


Bicalutamide 


4-(2,3-Dihydro-l  H-perimidin-2- 
yl)benzene-l  ,2-diol 


2-[(2-Phenoxyethyl)sulfanyl]-lH- 

benzimidazole 


Tamoxifen 


Fulvestrant 


CGP  52608 


SR1078 


AC1M6FQE 


T0901317 


-  AR:  LBS  antagonist  [31] 

-  DB08899  (DrugBank) 

-  Approved,  prostate  cancer 


-  AR:  LBS  antagonist  [32] 

-  DB01128  (DrugBank) 

-  Approved,  prostate  cancer 


-  AR:  AF-2  antagonist  [35] 

-  CID  761631  (PubChem) 

-  Experimental,  prostate  cancer 

-  IC50  =  18.64  pM  (eGFP  transcriptional  assay) 


-  AR:  BF-3  antagonist  [37] 

-  CID  1944110  (PubChem) 

-  Experimental,  prostate  cancer 

-  IC50  =  4.2  pM  (eGFP  transcriptional  assay) 


-  ERs:  [33]  LBS,  AF-2  antagonist 

-  DB00675  (DrugBank) 

-  Approved,  breast  cancer 


-  ERs:  LBS  antagonist  [34] 

-  DB00947  (DrugBank) 

-  Approved,  breast  cancer 


-  RORa  agonist  [69] 

-  CID  6509863 (PubChem) 

-  Experimental 


-  RORa  and  RORy  agonist  [71],  [72] 

-  CID  17980288  (PubChem) 

-  Experimental 


-  RORa  antagonist  (PubChem  BioAssay  AID:  610) 

-  CID  2347892  (PubChem) 

-  Experimental 

-  IC50  <  5.03  nM  (luciferase  reporter  assay) 


-  LXRs:  LBS  agonist  [93] 

-  CID  447912  (PubChem) 

-  Experimental 

-  EC50  =10  nM  (fluorescence  polarization  assay),  LXRa  [237] 


( continued  on  next  page) 
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Table  3  ( continued ) 
Compound  name 
GW3965 


Pyridaben 


XCT790 


DY131 


Diethylstilbestrol 


Genistein 


SRI 6388 


DIM-C-pPhOCH3  (C-DIM) 


Cytosporone  B 


Chemical  structure 


Annotation  (protein  target/accession  number/status/activity) 


-  LXRs  agonist  [88] 

-  CID  447905  (PubChem) 

-  Experimental 

-  EC50  =  40  nM  (FRET  based  cell-free  ligand  sensing  assay),  LXRp 
(PubChem  BioAssay  AID:  437600) 


-  COUP-TFII  antagonist  (PubChem  BioAssay  AID:  720548) 

-  CID  91754  (PubChem) 

-  Experimental 

-  IC50  =  61.93  nM  (luciferase  reporter  assay) 


-  ERRot  inverse  agonist  [113] 

-  CID  6918788  (PubChem) 

-  Experimental 

-  IC50  =  0.37  pM  (GAL4-ERRa  cell-based  transfection  assay)  [130] 


-  ERRp/y  agonist  [114,115] 

-  CID  5497124  (PubChem) 

-  Experimental 

-  EC50  =  0.13  pM  (FRET  assay),  ERRy  (PubChem  BioAssay  AID: 
246440) 

-  ERRy  antagonist  [128] 

-  DB00255  (DrugBank) 

-  Approved,  menopausal/postmenopausal  disorders 


-  ERRs  agonist  [29] 

-  CID  5280961  (PubChem) 

-  Natural  compound 


-  ERRot  antagonist  [131] 

-  CID  54612678  (PubChem) 

-  Experimental/Clinical  trial 


-  NGFIB  agonist/antagonist  [133] 

-  CID  11371604  (PubChem) 

-  Experimental 


-  NGFIB  agonist  [142] 

-  CID  10687292  (PubChem) 

-  Natural  compound 

-  EC50  =  0.115  nM  (luciferase  reporter  assay),  PubChem,  BioAssay 
AID:  500599 


6-Mercaptopurine 


-  NGFIB-like  agonist  [148] 

-  DB01033  (DrugBank) 

-  Approved,  antineoplastic  agent 
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Table  3  ( continued ) 


Compound  name 

Chemical  structure 

Annotation  (protein  target/accession  number/status/activity) 

Thiocolciran 

Q 

njT\-(jrV 

-  DAX1  antagonist  (PubChem  BioAssay  AID:  687017) 

-  CID  426728  (PubChem) 

-  Experimental 

-  IC50  =  2  1  2.03  nM  (luciferase  reporter  assay) 

Accession  numbers  are  provided  for  either  DrugBank  [6]  or  PubChem  [224],  and  chemical  structures  were  visualized  by  using  MarvinSketch. 


232,340  newly  diagnosed  cases  and  39,620  deaths  in  the  United 
States  [50].  Based  on  the  gene  expression  profiles,  breast  cancers 
can  be  divided  into  five  subtypes  that  vary  in  their  treatment 
options  and  survival  outcomes  [55-58].  ERa  positive  (ERa+)  and 
progesterone  receptor  positive  (PR+)  tumors  account  for  approxi¬ 
mately  70%  of  all  cases  [59,60  .  Additionally,  based  on  the  expres¬ 
sion  or  gene  amplification  of  the  ErbB-2/HER2  gene,  these  ERa+/ 
PR+  tumors  can  be  further  classified  into  HER2-  and  HER2+ 
subtypes.  For  both  early  and  advanced  ERa+  breast  cancer,  the 
standard  treatment  is  to  block  the  production  of  estrogens  (e.g. 
aromatase  inhibitors)  or  estrogen  binding  to  the  receptor  (e.g. 
Tamoxifen).  Adjuvant  Trastuzumab  therapy,  which  specifically  tar¬ 
gets  the  HER2  receptor,  is  used  to  treat  ERa-/PR-/HER2+  tumors 
and  has  resulted  in  positive  survival  outcomes.  In  contrast,  the 
ERa-/PR-/HER2-  subtype  or  triple  negative  breast  cancer,  which 
accounts  for  10-17%  of  all  breast  cancer  cases,  and  the  normal 
breast-like  (basal-like)  cancer  subtype,  which  accounts  for  15%  of 
the  cases,  are  both  without  approved  targeted  therapies  [61]. 
Unfortunately,  drug  resistance  to  even  the  most  effective  therapeu¬ 
tics  is  a  major  issue  affecting  breast  cancer  patients  today.  For 
example,  although  initially  effective,  one  third  of  patients  treated 
with  Tamoxifen  develop  recurrent  diseases  within  15  years  due 
to  the  progression  of  surviving  tumor  cells  to  a  hormone-resistant 
state  [62].  For  these  reasons,  alternative  drug  treatments,  which 
can  be  used  either  alone  or  in  combination  with  existing  therapies, 
are  needed. 

RAR-related  orphan  receptors  (RORs),  NR1F 

The  retinoic  acid-receptor-related  orphan  receptor  (ROR)  family 
consists  of  three  members:  RORa  (NR1F1),  RORJ3  (NR1F2),  and 
RORy  (NR1F3).  RORs  have  been  shown  to  regulate  lipid  and  steroid 
metabolism,  immune  response  and  circadian  rhythms.  RORa  is 
widely  expressed  in  many  tissues.  RORp  is  expressed  in  the  central 
nervous  system  in  particular  it  is  localized  to  regions  that  process 
sensory  information,  such  as  the  retina  and  pineal  gland  [63,64]. 
RORy  is  also  broadly  expressed  but  in  particular  it  is  expressed 
in  immune  tissues  such  as  the  thymus,  and  is  a  significant  factor 
in  the  development  of  T  helper  17  cells  65].  A  number  of  putative 
natural  and  synthetic  ligands  have  been  identified  that  target 
RORs;  most  notably  cholesterol,  cholesterol  metabolites  and  oxys- 
terols  that  have  been  shown  to  bind  to  the  LBD  of  RORs  and  induce 
conformational  changes  that  affect  coregulator  interactions 
[66,67].  While  recent  research  has  focused  on  the  therapeutic  value 
of  inhibiting  RORs  for  treating  metabolic  and  autoimmune  dis¬ 
eases,  analyses  of  published  gene  expression  data  indicate  that 
RORa  is  downregulated  in  multiple  cancers  including  prostate, 
breast  and  pancreatic  cancer,  and  suggest  that  RORa  may  act  as  a 
tumor  suppressor  [66,68-70]. 

In  prostate  cancer,  the  role  of  RORa  as  a  tumor  suppressor  was 
demonstrated  in  vitro  and  in  vivo  using  the  thiazolidinedione- 
derivative  CGP  52608  (Pubchem  CID:  6509863,  Table  3)  [69,70]. 
The  activation  of  RORa  by  CGP  52608  resulted  in  reduced  cell 


proliferation,  migration  and  invasiveness  in  vitro,  as  well  as 
reduced  growth  of  xenograft  tumors  in  nude  mice  using  the  andro¬ 
gen-independent  prostate  cancer  cell  line,  DU145  [69,70].  In  addi¬ 
tion,  these  studies  demonstrated  that  activation  of  RORa  leads  to 
an  increased  expression  of  the  cyclin-dependent  kinase  inhibitor 
p21,  but  decreases  the  expression  of  cyclin  A.  These  findings,  com¬ 
bined  with  the  observations  that  there  are  response  elements  for 
RORa  at  promoter  regions  upstream  of  several  cell  cycle-related 
genes  [69],  suggests  that  RORa  plays  an  inhibitory  role  in  cell  cycle 
regulation.  In  addition,  it  has  recently  been  shown  that  activation 
of  RORa,  as  well  as  RORy,  by  the  agonist  SRI 078  (Pubchem  CID: 
17980288,  Table  3)  can  stabilize  p53  protein  a  well-known  tumor 
suppressor;  an  observation  that  further  substantiates  a  tumor  sup¬ 
pressing  role  of  RORa  [10,71,72].  However,  further  investigation  is 
required  to  determine  if  activation  of  reduced  protein  levels  of 
RORa  in  cancer  would  be  enough  to  initiate  adequate  tumour 
suppression. 

In  breast  cancer,  several  lines  of  evidence  also  support  the  role 
of  RORa  as  a  tumor  suppressor.  RORa  expression  is  downregulated 
in  breast  cancer  tumor  cells  compared  to  normal  mammary  tissue 
[73,74].  Inhibition  of  RORa  expression  results  in  the  disruption  of 
polarized  acinar  cells  of  mammary  gland,  and  this  disruption  is 
known  to  be  an  early  event  in  breast  cancer  development  [75  .  In 
contrast,  overexpression  of  RORa  in  MCF12F  breast  cancer  cells 
resulted  in  a  reduced  rate  of  cell  proliferation  [73].  RORa  has  been 
shown  to  act  upon  several  pathways  that  are  known  to  be  involved 
in  breast  cancer  development  and  progression.  RORa  binds  to  (3- 
catenin  that  suppresses  Wnt/pcatenin  gene  transcription,  which 
is  known  to  be  an  important  pathway  in  cell  adhesion  and  cell 
growth  [76,77].  In  addition,  RORa  has  been  shown  to  have  a  nega¬ 
tive  effect  on  the  NF-kP  pathway,  which  is  known  to  regulate  genes 
important  in  cell  proliferation,  differentiation  and  apoptosis  in 
breast  cancer  [78-81].  Additionally,  RORa  plays  a  role  in  steroido¬ 
genesis.  In  RORa  knockout  mice,  the  regulation  of  hydroxysteroid 
dehydrogenases  &  sulfotransferases  were  disrupted  [82].  Further¬ 
more,  RORa  has  been  shown  to  enhance  the  expression  of  aroma¬ 
tase  in  ERa+  T47D  and  MCF7  breast  cancer  cells  through  its  binding 
to  the  aromatase  promoter  [83]. 

In  addition  to  the  discovery  of  compounds  CGP  52608  and 
SRI  078,  there  have  been  efforts  to  screen  for  RORa  inhibitors  at 
the  Scripps  Research  Institute  Molecular  Screening  Center  (Pub¬ 
chem  BioAssay  AID:  610).  They  reported  on  several  inhibitory 
compounds  with  IC50  values  of  <5  nM  in  a  transcriptional  activity 
assay  (AC1M6FQE,  Table  3).  There  are  currently  over  590  experi¬ 
mental  compounds  in  the  BindingDB  that  target  the  RORa, 
although  none  of  them  have  been  clinically  approved  (Fig.  1)  [7]. 
Based  on  the  existing  knowledge  on  natural  and  synthetic  ROR 
ligands  (reviewed  by  Kojetin  and  Burris  [64])  and  considering  the 
availability  of  a  crystal  structure  of  the  ligand  binding  domain  of 
RORa  in  complex  with  cholesterol  (PDB  ID:  1N83)  [84],  structural 
analyses  of  the  LBS  and  virtual  compound  screening  could  be  used 
to  discover  novel  RORa  agonists  or  antagonists  and  to  improve  the 
drug-like  properties  of  existing  candidates. 
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Liver  X  receptors  (LXRs),  NR!  FI 

As  the  name  suggests,  liver  X  receptors  (LXRs)  were  initially  iso¬ 
lated  from  liver,  but  subsequent  studies  have  shown  that  while  the 
expression  of  LXRoc  is  more  restricted  to  certain  tissues,  such  as  the 
liver,  kidney,  and  adipose,  its  p  isoform  is  ubiquitously  expressed  in 
almost  all  tissue  types  (Table  1)  [85].  The  main  functional  role  of 
LXRs  is  to  regulate  lipid  homeostasis,  and  activation  of  LXRs  leads 
to  the  reduction  of  intracellular  cholesterol  [86].  Since  cholesterol 
is  a  precursor  needed  for  androgen  synthesis,  activating  LXRs  can 
reduce  androgen  levels  in  prostate  cancer  cells  [86]. 

It  has  been  reported  that  a  synthetic  LXR  agonist  T0901317 
(Pubchem  CID:  447912,  Table  3),  which  mimics  the  structure  of 
the  endogenous  ligand  oxysterol,  was  able  to  decrease  cell  prolifer¬ 
ation  of  the  androgen-dependent  LNCaP  prostate  cancer  cell  line, 
as  well  as  suppress  the  growth  of  LNCaP  tumor  xenografts  in  nude 
mice  [85,87].  This  anti-proliferative  effect  of  the  synthetic  LXR  ago¬ 
nist  was  deemed  to  be  a  result  of  increased  expression  of  cyclin- 
dependent  kinase  inhibitor  p27  and  a  decreased  cell  population 
in  S-phase  [87].  In  addition,  a  recent  study  has  reported  that  acti¬ 
vation  of  LXR  by  agonist  GW3965  (Pubchem  CID:  447905,  Table  3) 
reduced  cellular  proliferation  and  migration  in  LNCaP  cells  and 
reduced  tumor  volumes  in  LNCaP-xenografted  nude  mice,  possibly 
through  regulation  of  the  suppressor  of  cytokine  signaling  3 
(SOCS3)  pathway  [88]. 

It  should  be  noted,  that  LXRs  are  both  anti-proliferative  and  lip- 
ogenic  factors  in  breast  cancer.  By  using  a  microarray  gene  analysis 
of  four  breast  cancer  cell  lines,  it  was  shown  that  LXR  ligands  dis¬ 
rupt  cell  proliferation  by  down-regulating  genes  in  the  E2F  tran¬ 
scription  factor  family  that  are  critical  for  cell  cycle  progression 
[89  .  Additionally,  LXR  induces  the  expression  of  key  regulatory 
genes,  such  as  fatty  acid  synthase,  leading  to  increased  triglyceride 
production  within  MCF7  breast  cancer  cells  90].  Furthermore,  tar¬ 
geting  LXR  also  stimulates  the  efflux  of  cholesterol  within  the  cell 
which  leads  to  the  inhibition  of  cell  proliferation  and  stimulation  of 
apoptosis  in  MCF7  cells  [91,92].  Although  the  full  functionality  of 
LXRs  and  their  ligands  are  yet  to  be  fully  understood,  with  the 
available  evidence,  LXRs  could  be  a  promising  target  for  the  treat¬ 
ment  of  breast  cancer. 

LXRs  are  one  of  the  most  challenging  members  of  the  ONR  fam¬ 
ily  to  target  on  account  of  their  essential  role  in  hepatic  lipogene- 
sis.  As  a  consequence  there  is  an  additional  layer  of  complexity  to 
the  rational  drug  design  in  that  the  therapeutics  must  be  tissue 
selective  to  reduce  side  effects.  There  are  no  approved  substances 
reported  in  the  DrugBank  [6]  despite  the  fact  that  there  are  over 
700  experimental  LXR  ligands  deposited  in  the  BindingDB  [7].  A 
number  of  crystal  structures  of  the  LXR  a  and  (3  forms  have  been 
resolved,  including  LXRp  in  complex  with  an  agonist  T0901317 
(PDB  ID:  1UPW)  [93].  Both  T0901317  and  GW3965  are  agonists 
for  both  LXRoc  and  LXRp.  While  the  development  of  selective  acti¬ 
vators  for  LXRoc  and  LXRp  has  been  difficult  due  to  their  highly  sim¬ 
ilar  amino  acid  sequence  (74%),  several  recent  studies  have 
described  compounds  that  display  partial  selectivity  to  LXRp  [94- 
97].  Overall,  although  challenging,  the  abundance  of  existing  com¬ 
pounds  and  crystal  structures  should  facilitate  both  ligand-  and 
structure-based  rational  drug  design  efforts  against  LXRoc  and 
LXRp. 

Chicken  ovalbumin  upstream  promoter  transcription  factors  (COUP- 
TFs),  NR2F 

Currently,  chicken  ovalbumin  upstream  promoter  transcription 
factor  (COUP-TF)  I  and  II,  (reviewed  by  Lin  et  al.  [98])  have  no 
endogenous  ligands  identified.  COUP-TFI  and  II  are  known  to  have 
important  roles  in  embryonic  and  organ  development  [99].  COUP- 
TFs  bind  to  DNA  repeats  as  homo-  or  heterodimers  [100,101]. 


COUP-TFs  interact  with  the  same  ERoc  and  ERp  DNA  response  ele¬ 
ments,  the  estrogen  response  element  (ERE);  however,  COUP-TFs 
bind  to  half  EREs  whereas  ERs  preferentially  bind  to  the  consensus 
EREs  [102].  Klinge  et  al.  demonstrated  that  COUP-TFs  are  able  to 
directly  interact  with  and  disturb  ER-DNA  binding  and  as  a  conse¬ 
quence  regulate  the  transcriptional  activity  of  ERoc  [103].  Although 
the  full  relationship  with  ER  is  still  unclear  as  a  study  by  Metivier 
et  al.  demonstrated  that  the  formation  of  a  COUP-TF/ERoc  complex 
resulted  in  enhanced  ERoc  transcriptional  activity  [104]. 

In  the  prostate,  COUP-TFs  may  adopt  either  a  tumor-promoting 
or  tumor-suppressing  role,  depending  on  the  cellular  context  and 
surrounding  tumor  microenvironment  [99].  On  one  hand,  it  has 
been  shown  that  COUP-TFII  can  promote  prostate  tumor  growth 
by  inhibiting  the  TGF-p-induced  growth  barrier  and  that  overex¬ 
pression  of  COUP-TFII  can  promote  tumor  metastasis  when  cou¬ 
pled  with  PTEN  deletion  [105].  On  the  other  hand,  another  study 
reported  that  COUP-TFII  can  repress  AR  activation  by  binding  to 
both  the  DBD  and  LBD  of  the  AR  and  that  overexpression  of 
COUP-TFII  reduces  proliferation  of  LNCaP  cells  [106].  In  addition, 
recent  genome-wide  gene  expression  analysis  and  transcription 
factor  binding  experiments  revealed  that  COUP-TFI  is  a  direct  tar¬ 
get  gene  of  the  AR  and  that  its  expression  is  downregulated  upon 
AR  binding  at  the  promoter  region  107].  This  study  further 
reported  an  inverse  relationship  between  COUP-TFI  gene  expres¬ 
sion  and  AR  activation  in  a  normal  prostate  development  phase. 

While  the  exact  function  of  COUP-TFII  in  breast  cancer  remains 
unclear,  its  mRNA  expression  in  ERa+  (MCF7,  T47D,  and  ZR75-1) 
and  ERoc-  (MDA-MB-231)  breast  cancer  cells  was  shown  to  be  sig¬ 
nificantly  increased  in  comparison  to  normal  epithelial  cell  lines 
[108].  Interestingly,  COUP-TFII  levels  were  reduced  in  Tamoxifen- 
resistant  lines  derived  from  MCF7  cells,  and  re-expression  of 
COUP-TFII  restored  the  growth  inhibitory  effect  of  Tamoxifen 
[104,109].  Further  to  this  study,  Litchfield  et  al.  suggested  that 
reduced  expression  of  COUP-TFII  in  the  endocrine  resistant  cell  line 
LCC9  contributed  resistance  due  to  the  suppression  of  NFkB  activ¬ 
ity  [110]. 

To  date,  there  are  a  number  of  compounds  reported  for  COUP- 
TFII  (but  not  for  COUP-TFI)  and  no  approved  drugs  that  target 
the  COUP-TF  family.  There  are  231  experimental  compounds  with 
reported  binding  affinities  toward  COUP-TFII  in  the  BindingDB 
deposit  [7].  Among  these,  pyridaben  (PubChem  CID:  91754,  Table  3) 
has  an  IC50  of  62  nM  that  was  determined  in  a  transcriptional  assay 
by  the  Scripps  Research  Institute  Molecular  Screening  Center  (Pub¬ 
chem  BioAssay  AID:  720548).  The  crystal  structure  of  the  ligand 
binding  domain  of  COUP-TFII  is  available  (PDB  ID:  3CJW)  and 
shows  that  oe-helix  10  extends  into  the  LBD  and  that  the  AF-2  is 
folded  into  a  position  that  blocks  coregulator  interactions  resulting 
in  autorepression  of  the  receptor  [111].  However,  the  receptor  still 
contains  a  viable  ligand  binding  pocket  that  when  bound  with 
ligand  releases  the  autorepression  [111].  Due  to  a  high  sequence 
similarity  (96%)  at  LBD,  the  COUP-TFII  structure  can  be  used  as  a 
template  for  COUP-TFI  for  further  in  silico  screening.  In  particular, 
in  silico  screening  should  be  applied  to  discover  chemical  struc¬ 
tures  with  high  binding  affinity  at  the  LBS  of  COUP-TFs.  Such  com¬ 
pounds  may  act  as  either  agonists  or  antagonists  and  can  be  used 
as  chemical  probes  to  further  elucidate  the  functions  of  COUP- 
TFs  in  prostate  and  breast  cancer  cells. 

Estrogen  related  receptors  (ERRs),  NR3B 

The  estrogen-related  receptors  (ERRs),  ERRoe,  ERRp  and  ERRy, 
were  so  named  for  their  sequence  and  structural  similarity  to 
estrogen  receptors  (ERs)  [112  .  While  initially  thought  to  share  a 
common  biological  function  with  the  ERs,  ERRs  do  not  bind  to 
estrogens  or  endogenous  ER  ligands.  There  are  no  known  endoge¬ 
nous  ligands  for  ERRs,  and  the  receptors  are  constitutively  active 


M.  Roshan-Moniri  et  ail  Cancer  Treatment  Reviews  40  (2014)  1 137-1152 


1147 


due  to  the  structure  of  the  LBD  that  leaves  the  AF-2  site  exposed 
and  capable  of  binding  coregulators  in  the  absence  of  ligand  bind¬ 
ing.  Thus,  the  transcriptional  properties  of  the  ERRs  are  mediated 
by  either:  coactivators,  such  as  peroxisome-proliferator  activated 
receptor  y  coactivator- lot  (PGC-loe)  and  the  pi 60  family  of  steroid 
receptor  coactivators  (SRC-1,  SRC-2,  SRC-3),  or  corepressors  such 
as  nuclear  receptor  interacting  protein  140  (RIP140).  The  active 
ERR  is  able  to  bind  to  DNA  at  the  5'-TNAAGGTCA-3'  sequence 
(where  N  is  any  nucleotide)  or  to  the  classical  ER  response  element 
(ERE)  S'-AGGTCANNNTGACCT-S'.  ERRs  regulate  genes  involved  in 
the  cellular  energy  metabolism  pathways,  including:  glycolysis, 
tricarboxylic  acid  cycle  and  oxidative  phosphorylation  [112].  Thus, 
ERRoe  is  predominantly  expressed  in  tissues  that  have  high  energy 
demand,  including  skeletal  muscle,  heart,  kidney,  liver  and  adi¬ 
pose.  Whereas,  ERRp  is  expressed  during  early  embryonic  develop¬ 
ment  and  ERRy  is  expressed  in  the  central  nervous  system  and 
spinal  cord. 

Of  the  ERRs,  ERRoe  is  known  to  promote  cell  proliferation  in 
multiple  cancer  cells  lines,  including  prostate  and  breast  cancer 
cells.  Furthermore,  inhibition  of  ERRoe  with  the  inverse  agonist 
(i.e.  binds  to  the  same  receptor  protein  but  induces  an  opposite 
response  to  an  agonist)  XCT790  (Pubchem  CID:  6918788,  Table  3) 
reduces  cell  proliferation  [113].  However,  overexpression  of  either 
ERRp  or  ERRy  resulted  in  the  suppression  of  cell  proliferation  in 
both  androgen-sensitive  and  androgen-insensitive  prostate  cancer 
cells,  and  this  growth  inhibition  could  be  further  enhanced  by  the 
ERRp/y  agonist,  DY131  (Pubchem  CID:  5497124,  Table  3)  [114- 
116].  In  addition,  transcription  reporter  and  chromatin  immuno- 
precipitation  assays  demonstrated  that  ERRp  directly  transacti- 
vates  a  promoter  upstream  of  the  cyclin-dependent  kinase 
inhibitor  p21  gene,  which  inhibits  cell  cycle  progression  [115]. 
While  the  mechanism  through  which  the  three  ERR  isoforms 
induce  their  effects  on  cell  proliferation  is  undetermined,  these 
findings  suggest  that  in  prostate  cancer,  ERRoe  acts  as  a  tumor  pro¬ 
moter,  whereas  ERRp  and  ERRy  act  as  tumor  suppressors. 

Due  to  their  structural  similarity  to  ERs,  ERRs  (in  particular 
ERRoe  and  ERRy)  are  most  recognized  for  their  potential  function 
in  breast  cancer.  Members  of  the  ERR  family  can  bind  to  EREs 
and  affect  gene  expression  of  ER-regulated  genes  and  may  contrib¬ 
ute  to  breast  cancer  development  or  progression.  Interestingly, 
Stein  et  al.  showed  that  the  gene  transcription  products  of  ERRoe 
had  little  overlap  with  those  of  ER  [117].  Furthermore,  the  roles 
of  both  ERRoe  and  ERRy  are  distinct  between  ERoe+  and  ERoe-  breast 
cancer.  ERRoe  can  act  as  an  activator  or  repressor  of  transcription  in 
ERoe-  or  ERoe+  cells,  respectively  [118].  For  example,  ERRoe  can 
increase  the  expression  of  the  estrogen-regulated  trefoil  factor  1 
(TFF1)  gene  in  ERoe-  breast  tumors  [119,120].  Similarly,  in  ERoe- 
breast  cancer  samples,  overexpression  of  ERRoe  induced  vascular 
endothelial  growth  factor  (VEGF)  [117,121].  In  comparison,  ERRy 
is  predominantly  associated  with  ERoe+/PR+  breast  tumors  in 
which  it  promotes  mesenchymal-to-epithelial  transition  (MET) 
by  stimulating  E-cadherin  expression,  which  results  in  an  inhibi¬ 
tion  of  tumor  growth  [122,123].  Furthermore,  ERRy  expression 
may  be  a  marker  of  tamoxifen  resistance  and  it  has  been  shown 
that  an  AAAG  tetranucleotide  repeat  polymorphism  in  untrans¬ 
lated  region  of  the  ERRy  gene  is  associated  with  breast  cancer  sus¬ 
ceptibility  [124,125].  While  the  exact  mechanism  of  action  of  both 
ERRoe  &  ERRy  in  breast  cancer  is  still  not  fully  understood,  current 
findings  support  a  potential  role  of  these  receptors  as  therapeutic 
targets  in  breast  cancer. 

There  are  relatively  few  compounds  reported  to  interact  with 
ERRs  (reviewed  Lu  et  al.  [126]).  There  are  several  crystal  structures 
resolved  for  the  ligand  binding  domains  of  ERRa  and  ERRy.  No 
crystal  structure  is  available  for  the  LBD  of  ERRp,  but  the  high  pro¬ 
tein  sequence  identity  (80%)  with  ERRy  makes  homology  modeling 
feasible.  However,  the  small  size  of  the  ligand  binding  pocket  has 


made  it  difficult  to  directly  target  ERRs  for  pharmacologic  activa¬ 
tion  127].  Hence,  the  BindingDB  catalogues  only  91,  31,  and  29 
compounds  that  can  bind  to  ERRa,  p  and  y,  respectively  [7].  The 
FDA-approved  drugs  diethylstilbestrol  (DB00255,  Table  3)  and 
tamoxifen  have  been  shown  to  interact  with  ERRs  and  act  as  antag¬ 
onists  [128,129].  Whereas,  XCT790  (Pubchem  CID:  6918788, 
Table  3)  binds  selectively  to  the  LBD  of  ERRa  and  has  an  IC50  of 
0.37  pM  [130].  Methods  that  combined  structure-based  virtual 
screening  and  biological  functional  assays  have  identified  several 
isoflavones  including  Genistein  that  bind  at  the  LBS  of  the  ERRs 
[29].  SRI 6388  (Pubchem  CID:  54612678,  Table  3),  a  novel  steroidal 
antiestrogen,  disrupts  the  interaction  between  ERRa  and  its  coacti¬ 
vator  peroxisome-proliferator  activated  receptor  y  coactivator- la 
to  inhibit  ERRa  activity.  SRI  6388  has  an  IC50  of  0.2  pM  on  malig¬ 
nant  cell  growth  in  xenograft  models  and  has  successfully  com¬ 
pleted  phase  I  and  phase  II  clinical  trials  [131]. 

Nerve  growth  factor  tB-like  ( NGFIB ,  NURR 1 ,  NOR1),  NR4A 

The  nerve  growth  factor  IB-like  family  is  comprised  of  three 
members:  NGFIB  (NUR77  or  NR4A1 ),  NURR1  (NR4A2),  and  NOR1 
(NR4A3).  These  nuclear  receptors  are  unique  due  to  the  fact  that 
they  lack  the  LBS  within  their  LBD  [45].  Members  of  this  family 
regulate  a  wide  range  of  biological  processes  including  metabo¬ 
lism,  inflammation,  vascular  function,  steroidogenesis,  and  func¬ 
tions  of  the  central  nervous  system  [132,133].  NIGFB  is  known  to 
interact  with  specific  response  elements  5'-AAAGGTCA-3'  in  the 
form  of  a  monomer  or  homodimer  [134,135].  The  physiological 
role  of  NGFIB  is  not  yet  fully  understood,  but  several  studies  indi¬ 
cate  it  has  a  role  in  cancer  cell  death  and  survival  [136-139].  Inter¬ 
estingly,  NGFIB  is  found  to  have  pro-survival  benefits  while  in  the 
nucleus,  but  promote  apoptosis  through  Bcl-2/cytochrome  c  path¬ 
way  when  translocated  into  mitochondria  [140,141].  NGFIB  has 
been  previously  shown  to  exhibit  pro-oncogenic  activity  in  several 
cancers,  including  cervical,  lung,  colon  and  pancreatic  cancers.  A 
number  of  compounds,  including  l,l-bis(3/-indolyl-l-(p-substi- 
tuted  phenyl)methanes  (C-DIM;  Pubchem  CID:  54612678,  Table  3), 
which  can  activate  or  deactivate  NGFIB,  have  been  reviewed  [133]. 
In  addition,  cytosporone  B  has  been  previously  reported  as  a  natu¬ 
rally  occurring  agonist  for  NGFIB  and  was  capable  of  promoting 
apoptosis  and  reducing  xenograft  tumor  growth  by  inducing  NGFIB 
expression  and  translocation  to  mitochondria  [142]. 

With  regard  to  prostate  cancer,  a  recently  reported  immunohis- 
tochemistry  experiment  revealed  that  expression  of  NURR1  pro¬ 
tein  is  higher  in  prostate  cancer  tissues  compared  to  benign  ones 
[143].  Furthermore,  this  NURR1  expression  was  positively  corre¬ 
lated  with  several  clinical  indicators,  including  tumor  stage  and 
Gleason  score.  Knockdown  of  NURR1  protein  with  siRNA  resulted 
in  reduced  proliferation,  migration  and  invasion  as  well  as  induced 
apoptosis  in  prostate  cancer  cells.  Hence,  the  study  suggested  that 
NURR1  may  be  used  as  a  biomarker  and  therapeutic  target  in  pros¬ 
tate  cancer  [143]. 

In  breast  cancer,  various  diindolylmethanes  (DIMs)  have  dem¬ 
onstrated  agonistic  effects  toward  the  NGFIB  [144].  Chinthalapalli 
et  al.  demonstrated  that  in  the  presence  of  DIM,  MCF7  breast  can¬ 
cer  cells  that  express  NGFIB  have  a  reduction  in  their  proliferation 
capacity  and  increased  apoptosis  [145].  In  addition,  NGFIB  is 
known  to  regulate  and  affect  transcriptional  responses  in  other 
ONRs  and  NRs,  such  as  COUP-TFs  and  retinoid  X  receptors  (RXRs). 
Furthermore,  it  has  been  demonstrated  that  NGFIB  is  capable  of 
producing  a  heterodimer  with  RXR  that  induces  apoptosis  and 
growth  inhibition  in  MDA-MB-231  breast  cancer  cells  through  acti¬ 
vation  of  RARa  and  subsequent  pathways  [146]. 

There  is  not  yet  an  approved  drug  targeting  any  member  of  the 
NR4A  subfamily.  The  recent  research  efforts  dedicated  toward  the 
NGFIB  have  generated  281  experimental  compounds  targeting  this 
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protein.  There  are  currently  two  crystal  structures  available  for 
NGFIB  (PDB  ID:  2QW4  and  3V3E  [147]),  which  should  further  aid 
in  the  identification  or  modification  of  new  or  existing  compounds. 
In  contrast,  there  are  only  12  compounds  which  have  been  shown 
to  interact  with  NURR1 .  The  crystal  structure  of  the  NURR1  LBD 
(PDB  ID:  10VL  [45])  is  observed  to  be  filled  by  hydrophobic  substit¬ 
uents,  which  suggests  that  it  may  be  a  true  ONR  without  a  ligand 
[45].  Interestingly,  6-mercaptopurine  and  a  few  structural  ana¬ 
logues  have  been  shown  to  activate  the  NR4A  subfamily.  The 
action  of  6-mercaptopurine  occurs  through  the  NTD  in  an  AF-1 
dependent  manner  [148].  The  NTD  of  NRs  is  traditionally  thought 
of  as  being  structurally  unstable  and  thus  has  not  been  the  focus 
of  many  drug  targeting  approaches.  However,  given  the  identifica¬ 
tion  of  NTD  directed  compounds,  further  work  in  this  direction 
may  yield  more  promising  results  that  could  be  applied  to  other 
ONR  and  NR  targets. 

DAX1,  NR0B1 

While  the  physiological  significance  of  DAX1  in  females  is 
unclear,  knockout  mouse  studies  have  demonstrated  that  DAX1 
is  essential  for  the  maintenance  of  male  reproductive  functions 
[149].  Furthermore,  DAX1  is  uniquely  and  highly  expressed  in  ste¬ 
roidogenic  tissues  [150].  For  example,  DAX1  is  known  to  be  a  neg¬ 
ative  regulator  of  aromatase  mRNA  expression  in  ovarian  cells  that 
are  the  sites  of  androgen  and  estrogen  production  [151].  Addition¬ 
ally,  several  studies  suggest  that  the  primary  function  of  DAX1  is  as 
a  coregulator  protein,  since  it  is  capable  of  affecting  the  activity  of 
several  receptors,  including  an  ONR,  steroidogenic  factor-1  (SF-1), 
as  well  as,  ER,  PR  and  AR  [152]. 

While  not  detected  in  normal  prostate  tissue,  DAX1  is  expressed 
in  prostate  cancer  epithelium,  and  its  expression  is  inversely  corre¬ 
lated  with  Gleason  score  [153].  DAX1  has  been  shown  to  directly 
interact  with  the  AR,  inhibiting  its  transcriptional  activity  through 
disruption  of  the  AR  N-  and  C-terminal  interaction  [154].  DAX1, 
which  is  able  to  cycle  between  the  cytoplasm  and  the  nucleus,  is 
also  capable  of  tethering  the  AR  in  the  cytoplasm  or  nucleus  of 
the  cell,  further  restricting  AR  transcriptional  activity  [154]. 

A  recent  publication  by  Conde  et  al.  established  that  while 
DAX1  is  expressed  in  the  benign  form  of  breast  carcinoma  (in  situ), 
it  is  more  significantly  expressed  in  invasive  breast  cancer  [155]. 
There  was  also  a  report  of  a  direct  correlation  between  lymph  node 
metastasis  and  the  expression  of  DAX1  protein.  In  addition,  it  was 
observed  that  DAX1  has  also  been  shown  to  directly  interact  with 
AR,  ERa  and  ER(3  in  breast  cells,  but  has  only  been  shown  to  influ¬ 
ence  the  intracellular  positioning  of  AR  and  ER-p,  and  not  ER-a 
[155,156].  Thus,  DAX1  may  be  a  novel  coregulator  of  nuclear 
receptor  action  in  both  prostate  and  breast  tissues. 

DAX1  is  unique  among  other  members  of  the  NR  family  in  that 
it  lacks  a  classical  DBD  domain  and  it  has  an  N-terminal  domain 
that  contains  three  repeats  of  a  unique  70  amino  acid  motif 
[157].  Additionally,  DAX1  contains  a  potent  transcriptional  repres¬ 
sor  domain  in  its  C-terminus  [158,159].  Although  there  is  not  yet  a 
protein  structure  available  for  human  DAX1 ,  the  LBD  of  the  mouse 
homolog  has  been  resolved  (PDB  ID:  3F5C  [160])  and  established 
that  the  C-terminal  domain  does  not  contain  a  ligand  binding  site. 
The  structure  of  the  human  domain  shares  68%  identity  with  the 
mouse  so  it  is  possible  that  the  human  DAX1  LBD  surface  may  con¬ 
tain  a  pocket  that  could  be  targeted  by  small  molecules.  Addition¬ 
ally,  drugs  may  also  be  developed  that  target  alternative  surfaces  of 
DAX1,  such  as  areas  involved  in  coregulator  interactions  or  the  N- 
terminal  domain  because  of  its  unique  sequence  repeat.  Currently, 
there  are  over  200  compounds  indicated  in  BindingDB  to  bind  to 
DAX1  [7].  The  Scripps  Research  Institute  Molecular  Screening  Cen¬ 
ter  (PubChem  BioAssay  AID:  687017)  reported  that  compound  thi- 
ocolciran  (Pubchem  CID:  426728,  Table  3)  demonstrates  an  IC50  of 


212  nM  in  a  luminescence-based  high-throughput  cell  screening 
assay.  Additionally,  there  are  two  FDA-approved  drugs  that  have 
been  shown  to  indirectly  target  DAX1  [6].  One  is  dexamethasone, 
a  potent  glucocorticoid,  and  the  other  is  Tretinoin,  a  naturally 
occurring  derivative  of  vitamin  A.  Thus,  DAX1  is  a  promising  target 
due  to  its  unique  structure  that  may  provide  an  alternative  strat¬ 
egy  for  drug  targeting. 

Conclusions 

In  recent  decades  nuclear  receptors  have  been  actively 
exploited  as  drug  targets  mainly  due  to  their  regulatory  roles  in 
numerous  cellular  processes  and  the  presence  of  “druggable” 
ligand  binding  domain.  Unlike  most  other  transcriptional  factors 
that  lack  a  well-defined  active  site,  the  function  of  NRs  can  be  mod¬ 
ulated  by  small  molecules  that  compete  with  endogenous  ligands 
to  act  as  either  antagonists  or  agonists.  Within  the  NR  family,  hor¬ 
mone  receptors  such  as  the  AR  and  ER  have  been  studied  inten¬ 
sively,  and  represent  established  targets  for  a  number  of 
approved  drugs  in  prostate  and  breast  cancers.  However  for  both 
cancers,  current  hormonal  therapies  often  lead  to  the  development 
of  drug  resistance;  thus,  targeting  alternative  proteins  and  path¬ 
ways  may  provide  a  way  to  treat  resistant  diseases.  Given  the 
recent  advances  in  computational  modeling  techniques  and  the 
availability  of  protein  structures,  ONRs  present  an  exciting  oppor¬ 
tunity  for  drug  development.  As  described  in  this  review,  with  the 
aid  of  in  silico  modeling  much  progress  has  been  made  on  the  iden¬ 
tification  of  natural  and  synthetic  ONR  ligands.  Furthermore  the 
recent  successes  targeting  alternative  protein-protein  interaction 
sites  such  as  the  AF-1,  AF-2,  and  BF-3  site  of  the  AR  highlight  the 
possibilities  for  targeting  these  sites  on  ONRs.  Ultimately,  identifi¬ 
cation  of  novel  small  molecules  that  bind  to  ONRs  will  provide  a 
better  understanding  of  their  biology,  structure  and  interactions. 
In  addition,  this  accumulated  knowledge  on  the  endogenous  func¬ 
tions  of  ONRs  will  provide  a  solid  foundation  for  future  drug  dis¬ 
covery  with  therapeutic  values  that  may  go  beyond  prostate  and 
breast  cancers. 
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